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ABSTRACT 

The Blanco Cosmology Survey (BCS) is a 60 night imaging survey of ~80 deg^ of the southern sky 
located in two fields: {ol^5)= (5 hr, —55°) and (23 hr, —55°). The survey was carried out between 
2005 and 2008 in griz bands with the Mosaic2 imager on the Blanco 4m telescope. The primary 
aim of the BCS survey is to provide the data required to optically confirm and measure photometric 
redshifts for Sunyaev-Zel'dovich effect selected galaxy clusters from the South Pole Telescope and the 
Atacama Cosmology Telescope. We process and calibrate the BCS data, carrying out PSF corrected 
model fitting photometry for all detected objects. The median lOcr galaxy (point source) depths over 
the survey in griz are approximately 23.3 (23.9), 23.4 (24.0), 23.0 (23.6) and 21.3 (22.1), respectively. 
The astrometric accuracy relative to the USNO-B survey is ~ 45 milli-arcsec. We calibrate our 
absolute photometry using the stellar locus in griz J bands, and thus our absolute photometric scale 
derives from 2MASS which has ~ 2% accuracy. The scatter of stars about the stellar locus indicates 
a systematics fioor in the relative stellar photometric scatter in griz that is ~1.9%, ~2.2%, ~2.7% 
and~2.7%, respectively. A simple cut in the AstrOmatic star-galaxy classifier spreadjnodel produces 
a star sample with good spatial uniformity. We use the resulting photometric catalogs to calibrate 
photometric redshifts for the survey and demonstrate scatter Sz/{1 -\- z) = 0.054 with an outlier 
fraction 77<5%toz~l. We highlight some selected science results to date and provide a full 
description of the released data products. 

Subject headings: cosmology:cosmic microwave background - cosmology: observations - galax- 
ies:clusters - Sunyaev-Zel'dovich Effect . Optical Surveys 



1. INTRODUCTION 



Since the discovery of cosmic acceleration at the end 



of the last millenium (Schmidt et al. 1998 
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et al. 1999"), understanding the underlying causes has 



remamed as one of the key mysteries in modern astro- 
physics. As the most massive collapsed structures in the 
universe, galaxy cluster populations and their evolution 
with redshift provide a powerful probe of, for example, 
the dark energy equation of state parameter as well as al- 
ternate gravity theo ri es wh i ch mimic cosmic acce l eration 



( Wang fc Steinhardt] [1998 
,et al.-2UUl|). Evolution olt 



Haiman et al. 2001 Holder 



le cluster abundance depends 
on a combination of the angular-diameter distance vs. 
redshift relation and the growth rate of density pertur- 
bations. This sensitivity enables one to constrain a range 
of cosmological parameters, inclu ding the ma tter density, 
the sum of the neutrino masses (Ichiki & Takada 2011| ), 
the present day amplitude of density fluctuations, and 
the presence of primord ial non-Gaussianity in the ini- 



tial d ensity fluctuations (Dalai et al.||2008 Cunha et al. 



2010). In addition, galaxy clusters provide a n ideal la bo- 



ratory to study galaxy evolution (e.g. Dressl er|1980 ). In- 



e.e 



Butcher 



teresting studies of the galaxy properties and their evo- 
lution within clusters include studies of the blue frac- 
tion and the ha l o occupation di s tribution 
Oemler"198 4| |Lin et al.||2 003; Tm 



. . Mohr " "2bjj4| ILm 

et al. 2006; Hansen et al.||2b09; Zenteno et al. 20TT| 

The first large scale attempt to identify and catalog 
galaxy clusters was by Abell in 1958. He discovered 



galaxy clusters by looking for over-densities of galax- 
ies in Palomar Observatory photographic plates within 
a radi us of about 2.1 Mpc around a given cluster po- 
sition (Abell 1958). Abell's catalogs contained about 



2 



Desai et al. 



4700 clusters ( |Abell et aL]|1989| ). However Abell's cata- 
log suffered from incompleteness and cont amination from 
projection effects as well as human bias ( Biviano||2008 ) . 
With the advent of CCD cameras one could apply ob- 
jective automated algorithms to look for galaxy clusters, 
and this has led to significant progress in cosmological as 
well astrophysical studies using galaxy clusters. 

In the last decade, many optical photometric surveys 
such as SDSS, CFHTLS, RCS covering contiguous re- 
gions of the sky have discovered several new galaxy clus- 
ters spanning a broad range of masses and redshifts. The 
CFHTLS- W dAdami et aH2010| ) has observed about 171 
deg^ in griz bands with 80 % completeness up to i band 
magnitude of 23. The RCS-2 ( [Gilbank et aL]|2QTT] ) sur- 
vey has covered approximately 1000 deg^ in grz bands 
with lOcr magnitude depths of around 23.55 in r-b and. 
The SDSS MaxBCG catalog ( [Koester et al.||2QQ7| has 
covered about 7,500 deg^ in ugriz bands and with 10 cr r 
band magnitude limits of about 22.35. The largest opti- 
cal galaxy cluster survey in terms of area is the Northern 
Sky Optical Cluster Survey HI which has imaged about 
11,40 deg^ up to a redshift of about 0.25 (Gal et al. 
|2009[ ). The deepest optical cluster sur vey to date is the 
lTFHTLS-D survey ( |Adami et al.|2010D , which reaches 80 
% completeness for i band magnitudes of 26 and has de- 
tected clusters up to redshift of 1.5. Two upcoming pho- 
tometric galaxy cluster surveys which will start around 
October 2012 include DES which will cover about 5,000 
deg^ in grizY bands wit h 10 a r-band limit ing magni- 
tudes of 24.8, and KIDS ( |de Jong et al.||2012D which will 
cover 1,500 deg^ in ugri bands with lOcr r-band limiting 
magnitude of 24.45. 

One can use such surveys for cos mological studies 
using galaxy clusters. For example, [Gladders et al. 



(|2007') showed that a large optical galaxy cluster sur- 
vey could constrain cosmologi cal parameters using th e 
self-calibration method (Majumdar fc Moh r|20Q3l 12004') . 
First cosmological con straints usi ng S DSS optical cata- 
logs are described in [Rozo et al.| (|2010| . 

Over the last decade there have been several mm- 
wave cluster stu dies in the Southern H emisphere , in- 
cluding ACBAR ("Reichar dt et al. I '20091, ACT (Fowler 
|004), APEX (Gonzalez et al.||2001) and SPT (Ruhl 
^t al. 2004). All these projects have attempted to carry 
out galaxy cluster surveys using the Sunyaev-Zel'dovich 
effect (SZE). The SZE is the distortion of cosmic mi- 
crowave background spectrum due to inverse Compton 
scattering of CM B photons by hot el ectrons in galaxy 
clusters (Sunyaev & Zerdovich"1972'), and it provides 
a promising way to discover galaxy clusters. Because 
the surface brightness of the SZE signature of a particu- 
lar cluster is independent of redshift, SZE survey cluster 
samples can in princi ple have sensitivity over a broad 
range of redshifts (Birkinshaw '1999( 'Carlstrom et al. 
[2002 ). However, to make use of SZE selected galaxy 
cluster samples, one needs a well-understood selection 
of galaxy clusters (sample contamination and complete- 
ness), cluster redshift estimates and a link between SZE 
signature and the cluster halo masses. It is important to 
note that redshift estimates cannot be obtained using SZ 
experiments alone, and so one needs dedicated optical 
surveys to follow up these galaxy clusters detected by SZ 
surveys. 



The Blanco Cosmology Survey (BCS) is an optical pho- 
tometric survey which was designed for this purpose and 
positioned to overlap the ACBAR, ACT, APEX and SPT 
surveys in the southern hemisphere. The goal of BCS is 
to enable cluster cosmology by providing the data to con- 
firm galaxy clusters from the above surveys and to mea- 
sure their photometric redshifts. This was done by sur- 
veying two patches totalling ~80 deg^ positioned so that 
they could be observed with good efficiency over the full 
night during the period October - December from Chile. 
The BCS observing strategy was chosen to obtain depths 
roughly two magnitudes deeper than SDSS, so that one 
could estimate photometric redshifts for L > galaxies 
out to a redshift z = 1. 

The outline of this paper is as follows: Sect. [2] describes 
the Blanco Cosmology Survey, including the camera, ob- 
serving strategy and site characteristics. In Sect. [3] we 
describe in detail the processing and calibration oi the 
dataset using the Dark Energy Survey Data Manage- 
ment System. In Sect. |3.31 we describe the photomet- 
ric characteristics of the BCS dataset and present single 
galaxy photometric redshifts that are tuned using fields 
containing large numbers of spectroscopic redshifts. In 
this paper all magnitudes refer to AB magnitudes. 

2. BCS SURVEY 

Blanco Cosmology Survey was a NOAO Large Sur- 
vey project (2005B-0043, PI: Joseph Mohr) which was 
awarded 60 nights between 2005 (starting from semester 
2005B) and 2008 on the Cerro Tololo Inter American 
Observatory (CTIO) Blanco 4 m telescope using the Mo- 
saic2 imager with griz bands. Because of shared nights 
with other programs, the data acquisition included 69 
nights, and the final processed dataset only consists of 
66 nights, because two nights were entirely clouded out 
and the pointing solution for one night (20071105) was 
wrong due to observer error. We now describe the Mo- 
saic2 imager on the Blanco telescope and then discuss 
the BCS observing strategy. 

2.1. Mosaic2 Imager 

The Mosaic2 imager is a prime focus camera on the 
Blanco 4m telescope that contains eight 2048x4096 CCD 
detectors. The 8 SITe 2K x 4K CCDs are read out in 
dual- amplifier mode, where different halves of each CCD 
are read out in parallel through separate amplifiers. The 
CCDs are read out through a single amplifier per chip si- 
multaneously to 8 controller inputs. Read noise is about 
6 — 8 electrons and readout time is about 110 seconds. 
The dark current rate is less than 1 electron/pixel/hour 
at 90 K. The resulting mosaic array is a square of about 
5 inches on an edge. The gaps between CCDs are kept 
to about 0.7 mm in the row direction and 0.5 mm in the 
column direction. Given the fast optics at the prime fo- 
cus on the Blanco, the pixels subtend 0.27^' on the sky. 
Total field of view is 36.8 arc-minute on a side for a to- 
tal solid angle per exposure of ~0.4 deg^. More details 
on the Mosaic2 imager can be found in the online CTIO 
documentation^^ . 

2.2. Field selection and multi-wavelength coverage 



http : //www . ctio . noao . edu/mosaic/manual/index . html 
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Fig. 1. — BCS survey footprint of coadded tiles in the 5 hr and 
23 hr fields. There are 104 tiles covering -^35 deg2 in the 23 hr 
field and 138 tiles covering ~45 deg^ in the 5 hr field for a total 
coverage of ~80 deg-^. The black vertically hatched boxes represent 
tiles which have passed our quality checks. The red horizontally 
hatched boxes represent tiles with some data quality problems that 
we have not corrected. 



The survey was divided into two fields to allow effi- 
cient use of the allocated nights between October and 
December. Both fields lie near 5 = —55° which allows for 
overlap with the SPT and other mm- wave surveys. One 
field is centered near a = 23.5 hr and the other is at a = 
5.5 hr. The 5 hr 30 min —52° patch consists of a 12 x 11 
array of Blanco pointings and the 23 hr —55° patch is 
a 10 X 10 array of pointings. The 5 hr field lies within 
the Boomerang field where the ACBAR experiment took 
data. The 23 hr field has been observed by the APEX, 
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Fig. 2. — Redshift evolution of a passively evolving L* galaxy 
along with target lOcr photometric BCS depths in each band. The 
exposure times in each band were tuned so that photometric depth 
meets or exceeds L* out to the redshift where the 4000Abreak 
shifts out of that band, but also limited to z = 1 due to the low 
sensitivity of the Mosaic2 camera in the z-band. 



ACT and SPT experiments. In addition to the large sci- 
ence fields, BCS also covers nine small fields that overlap 
large spectroscopic surveys so that photometric redshifts 
using BCS data can be trained and tested using a sample 
of over 5,000 galaxies with spectroscopic redshifts. BCS 
also surveyed standard star fields for photometric cali- 
bration. The coverage of BCS in 5 hr and 23 hr fields 
is shown in Figs [l] For convenience of data process- 
ing and building catalogs, we divide the survey region 
into 36' X 36' square regions called tiles. Each tile is a 
8192 X 8192 pixel portion of a tangent plane projection. 
These tiles are set on a grid of point separated by 34', 
allowing for approximately V overlaps of sky between 
neighboring tiles. The black vertical hatches in Fig. [l] 
indicate locations of tiles which passed various quality 
checks. The red horizontal hatches indicate locations of 
tiles which were observed and processed, but failed data 
quality checks. 

We also secured other multi-wavelength observations 
overlapping parts of the BCS fields. About 14 deg^ of 
the 23 hr BCS field was surveyed using XMM-Newton 
(known as XMM-BCS survey) and results from those ob - 



servations are reported elsewhere (Suhada et al. 2012). 
An ~12 deg^ region of the same field was also targeted 
in a Spitzer survey (S-BCS). More recently, the XMM- 
Newton survey has been expanded to 25 deg^, and the 
Spitzer survey has been expanded to 100 deg^. Most of 
the BCS region has been observed in the near - infrared as 
part of the the ESO VISTA survey program (Cioni et al. 
|2011|. 

2.3. Observing Strategy 
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The BCS observing strategy was designed to allow us 
to accurately measure cluster photometric redshifts out 
to redshift z = 1. Because the 4000 A break is redshift- 
ing out to 8000A by z = 1, obtaining reliable photo-z's 
for < z < 1 requires all four photometric bands r, 

z (i.e. one loses all clusters at z < 0.4 if you drop 
^-band and with z-band we can actually push beyond 
z = 1). The redshift at which the 4000 A break red- 
shifts beyond a particular band sets, crudely speaking, 
the maximum redshift for which that band is useful for 
cluster photo-z's; for griz this is z = 0.35, 0.7, 1.0 and 
1.4, respectively. Because the central wavelength of the 
g-band is about 4800A with a FWHM of 1537A we start 
losing sensitivity to very low redshift clusters, because it 
is not possible to straddle the 4000A break. Although 
detailed studies of the sensitivity of optical cluster de- 
tection at low redshifts have not been done, our ability 
to estimate unbiased red sequence redshifts for clusters 
is reduced below redshifts z ^ 0.1 . 

We calculate our photometric limits in each band by 
requiring that the depth allows us to probe at least to 
at that maximum redshift with lOcr photometry. We 
use a Bruzual and Chariot z = 3 single burs t model with 
passive evolution ( Bruzual fc Charlot||2003 ) to calculate 
the evolution of in the lour bands (see Fig. [2|. We 
select our z depth to probe to at z = 1 rather than 
at z = 1.4, because of the low sensitivity of the Mosaic2 
detectors in the z-band. The survey was designed to 
reach lOa photometric limits within a 2.2 arcsec aperture 
of ^ = 24.0, r = 23.9, i = 23.6 and z = 22.3. These limits 
assume an airmass of 1.3 and 0.9^' median seeing for all 
bands. Assuming bright time for z and i and dark time 
for g and r these limits require exposures of 250 s, 600 s, 
1400 s and 700 s in griz^ respectively. 

In all, we observed about 288 tiles spanning our survey 
fields. For each field we typically took 2 exposures in g 
of 125 sees each, 2 exposures in r of 300s each, 3 expo- 
sures in i of 450s each and three exposures in z of 235 
seconds each. A limitation of the Mosaic2 detector is a 
very low saturation of around 25,000 ADU for most of 
the detectors, and this forced us to take short exposures 
even through the readout time for each was quite high. 
Neighboring pointings have small overlaps, but multiple 
exposures were offset by approximately half the width of 
an amplifier to help us tie the survey together photomet- 
rically. Having two shifted exposures allows us to largely 
overcome the gaps in our survey left by spaces between 
neighboring chips. In addition to this primary survey 
tiling, we also constructed another layer of tilings which 
was designed to sit at the vertices of unique groups of four 
adjacent primary pointings. These tiles were observed 
using shorter exposures during poor seeing conditions on 
photometric nights. The 110 s readout of the Mosaic2 
camera makes the efficiency of short exposures low, and 
so in each band we have chosen the minimum number of 
exposures allowable given the sky brightness. The total 
exposure per tile is 3000 s and after including the readout 
time, the total observation time per science field is about 
4200 s, giving us an overall efficiency of about 70%. The 
dome ffats and bias frames were taken in the afternoon, 
and we did not take any twilight ffats. Over the course of 
the survey we acquired just over 3000 science exposures 
and an additional 455 photometric overlap exposures. 
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Fig. 3. — Sky brightness distributions for all four bands averaged 
on a per exposure basis during the BCS survey. Typically we ob- 
served in g and r during dark time and i and z during bright time. 
The brightness valuesare peaked at around 22.5, 21.5, 20.5, and 
18.75 mag/arcsec^ in griz bands, respectively. 

In addition to science exposures, on photometric nights 
we also observed photometric calibration fields as well as 
fields for calibrating our photo-z algorithms. These fields 
were CN0C2, DEEP, CFRS, CDFS, SSA22 and VVDS 
fields. For the photometric calibration fields we typically 
observed two or three fields during evening and morning 
twilight and a single field during the transition from the 
23 hr field to the 5 hr 30 min field. We observed in 
all four bands during these calibration exposures. The 
spectroscopic standard fields were observed to full science 
depth using the same strategy as for the full survey. 

2.4. Site Characteristics 

The BCS survey provides a sampling of the CTIO site 
characteristics over a 69 night period in the October to 
December timeframe over four observing seasons. Be- 
cause this is the same timeframe planned for DES ob- 
servations this provides an interesting glimpse into the 
expected site characteristics for DES. Given that the en- 
tire Mosaic2 camera and wide field corrector are being 
replaced by DECam and the new DECam wide field cor- 
rector (Soares-Santos et al. 2011), the seeing distribution 
for the DES data could be significantly improved relative 
to the BCS seeing distribution. 

The seeing distribution is shown in top panel of Fig. [Sj 
The seeing was obtained by running PSFEX software on 
all single-epoch images and using the FWHM_MEAN param- 
eter. The FWHM_MEAN is derived from elliptical Moffat 
fits to the non-parametric PSF models. These FWHMs 
include the pixel footprint. The modal seeing values in- 
tegrated over the survey are 1', 0.95^ 0.8^ 0.95' for 
griz bands, respectively. The median seeing values are 
1.07, 0.99, 0.95, and 0.95 arcseconds, while the upper and 
lower quartile seeing values are [0.96, 1.26], [0.89, 1.16], 
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Fig. 4. — The airmass distributions for BCS exposures, color 
coded by band and normalized by total number of exposures. The 
peak airmass values in griz bands are 1.144, 1.147, 1.138 and 1.141 
respectively. 



[0.84, 1.13], [0.83, 1.11] arcseconds respectively. 

The sky brightness is shown in Fig. [3j The sky bright- 
ness is calculated using ZP — 2.5 log 5, where ZP is the 
calculated zeropoint for that image and B is the sky 
brightness in ADU/arcsec^ . The sky brightness distri- 
butions in the griz bands have modal values of approx- 
imately 22.5, 21.5, 20.5 and 18.75 mag/arcsec^, respec- 
tively. Moreover, almost all i and z band data were taken 
with the moon up, while almost all g and r band data 
were taken with the moon set. The median values are 
22.3, 21.3, 20.3 and 18.7 mag/arcsec^, respectively. 

Given the division of the survey into a 23 hr and a 
5 hr field, it was possible to obtain most of the data at 
relatively low airmass. Fig.[4]shows the airmass distribu- 
tions for each band during primary survey observations. 
We often obtained photometric calibration field observa- 
tions over a broader range of airmasses, but we tried to 
restrict our primary survey observations to air masses of 
< 1.6. The median air mass in bands griz are 1.144, 
1.147, 1.138 and 1.141 respectively. 

3. DATA PROCESSING AND CALIBRATION 

The processing of BCS data is carried out using 
the automated Dark Energy Survey data management 
(DESDM) system which has been u nder development 



since Fall 2005 at University of Illinois (Ngeow et al. 2006 
|Mohr et al.]|2008D . DESDM wih be used to process, call 



brate and store data from the Dark Energy Survey once 
it begins operations in October 2012. Since 2005, the 
DESDM system has been validated through a series of 
data challenges with simulated DECam data, which en- 
abled us to improve various steps of the pipeline. The 
same automated pipeline was used to analyze BCS data. 
The only addition/change to the DESDM pipeline to an- 
alyze BCS data was in the crosstalk correction code, for 



which the routine had to be customized for the Mosaic2 
camera. Processing of the BCS data presented here has 
been carried out on National Teragrid resources at NCSA 
and LONI supercomputers together with dedicated work- 
stations needed for orchestrating the jobs and hosting 
the database. The middleware for the data reduction 
pipeline is designed using Condor batch processing sys- 
tem. Each night takes about 300 CPU hours for process- 
ing. 

We have processed BCS data multiple times in a pro- 
cess of discovery where we found problems with the data 
that required changes to our system. Scientific results 
from earlier rounds of processing of BCS data have al- 
ready appeared, including the optical confirmation of the 
first e ver SZE selected galaxy clusters (iStaniszewski et al. 
|2009 [) and the discovery of strong gravitational lensing 
arc using data from the first rou nd of processing in Spring 
2008 (Buckley-Geer et al.|2011[ ). Additional galaxy clus- 
ter science arising from subseque nt rounds of BCS pro- 
cessing have also been published ([High et ar||2010 Zen- 



|teno et al.|2QTT||Suhada et al.|2012|7" Currently our latest 
processing is being used for additional SZE cluster sci- 
ence within SPT, continued studies of the X-BCS region, 
and for the followup of the broader XMM-XXL survey 
over the 23 hr field. 

The BCS data were made public one year after their 
acquisition, as is standard policy at NO AO. This has 
enabled multiple independent teams to access the data 
and use it for their own scientific aims. The first three 
seasons of BCS data were processed using an i ndepen- 
dent pipeline developed at Rutgers University (Menan- 
teau et al | |2010bt i Menanteau fc Hughes||2Q09| . All tour 
seasons of BCS data have also been processed indepen- 
dently using NOAO pipeline as part of the current auto- 
mated proces sing program, an d with the PHOTPIPE anal- 
ysis pipeline ( Rest et ar||2005 ). 



3.1. Detrending 

In this section we describe in detail the key steps in- 
volved in the DESDM pipeline used for reduction of 
Mosaic2 data to convert raw data products to science 
ready catalogs and images. Data from every night 
are processed through a nightly processing or detrend- 
ing pipeline. Then data from different nights in the 
same part of the sky are combined using the co-addition 
pipeline. The detrending pipeline briefiy consists of 
crosstalk corrections, overscan, fiatfield, bias and illumi- 
nation correction, astrometric calibration and cataloging. 
We now describe in detail each step of the detrending 
pipeline. 

3.1.1. Crosstalk Corrections 

A common feature of multi-CCD cameras, such as Mo- 
saic2 imager, is crosstalk among the signals from oth- 
erwise independent amplifiers or CCDs. This leads to 
a CCD image containing not only the fiux distribution 
that it collected from the sky, but also a low amplitude 
version of the sky fiux distributions that appear in other 
CCDs. The crosstalk correction equation is described by: 



N 



i=i 



(1) 
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where li denotes the crosstalk corrected image pixel value 
in CCD, aij denote the cross-talk coefficients and /^w 
is the raw image pixel value. We used cross-talk coeffi- 
cients provided by NOAO through the survey. As part 
of the crosstalk-correction stage, the raw image (which 
contains 16 extensions) is split into one single-extension 
ffie per CCD. The processing and calibration of CCD mo- 
saics can proceed independently for each CCD after the 
crosstalk correction, and therefore we split the images 
to enable efficient staging of the data to the compute 
resources. 

3.1.2. Image Detrending 

Detrending is the process that removes the instrumen- 
tal signatures from the images. Detrending, in this con- 
text, includes overscan correction, bias subtraction, fiat 
fielding, pixel-scale correction, fringe and illumination 
correction. Both the overscan correction and bias cor- 
rection are required to remove the bias level present in 
the CCD and any residual, recurrent structure in the DC 
bias. Overscan correction is done for all raw science and 
calibration images. We subtract the median pixel value 
in the overscan region in each row for both the amplifiers 
in each CCD from the raw image pixel values after the 
crosstalk correction stage. 

The median bias frame is created using nightly bias 
frames taken during the late afternoon, and subtracted 
from the nightly data. The fiat field correction is typ- 
ically derived from dome fiats taken for each observing 
band. The input dome fiat images are overscan corrected, 
bias corrected, and then scaled to a common mode and 
then median combined. The resulting fiat field correc- 
tion is scaled by the inverse of the image mode, creat- 
ing a correction with mean value of about unity. For 
the bias correction and the fiat correction the variation 
among the input images is used to create an inverse vari- 
ance weight map that is stored as a second extension in 
the correction images. The creation of correction images 
also requires a bad pixel map, which is an image where 
pixels with poor response or with high dark current are 
masked and excluded from the images. These bad pixel 
maps are created initially using bias correction and fiat 
field correction images to identify the troublesome pixels. 

The bias and fiat field corrections are then applied to 
the science images to remove pixel to pixel sensitivity 
variations. These corrections are only applied to those 
science pixels that are not masked. In this process, 
each science image receives an associated inverse vari- 
ance weight map that encodes the Poisson noise levels 
and Gaussian propagated noise from each correction step 
on a per pixel basis. In addition, each science image has 
an associated bad pixel map (short integer) where a bit is 
assigned to each type of masking (i.e. pixels masked from 
the original bad pixel map, or masked due to saturation, 
cosmic ray, etc). In our data model the science image 
has three extensions: image, weight map, bad pixel map. 
Each measured fiux at the pixel level comes along with 
its statistical weight and a history of any masking that 
has been done on that pixel. 

For the Mosaic2 imager which has significant focal 
plane distortion, the pixel scale varies significantly over 
the field, leading to a significant trend in delivered pixel 
brightness as a function of position even with a fiat in- 
put sky. For such detectors flattening the sky introduces 



a photometric non-flatness to the focal plane. Typically 
this pixel scale variation is corrected during the process 
of remapping to a portion of a tangent plane, but in our 
case we prefer to do the single epoch cataloging on im- 
ages that do not suffer from correlated noise. Therefore, 
we apply a pixel scale correction to account for varia- 
tion of pixel response as a function of x and y position 
for each CCD. We first created master template images 
to determine photometric flatness corrections using as- 
trometrically reflned images from the Mosaic2 camera 
that we use to calculate the solid angle of each pixel. 
The correction image is then normalized by the median 
value, providing a flat fleld like correction image that 
can be used to bring all pixels to a uniform flux sensitiv- 
ity. To avoid reintroducing trends in the sky with this 
correction, we apply this correction only to the values 
of each pixel after subtraction of the modal sky value. 
Effectively, this correction scales only source flux while 
maintaining a flat sky. 

Illumination and fringe corrections are derived from 
fully processed science observations in a particular band. 
These can be from a single night or shared across nights. 
Usually if there was only one exposure from a given band 
in a night we use science observations from neighbouring 
nights to create the illumination and fringe correction 
images. Illumination corrections are done for all images, 
but fringe corrections on the Mosaic2 camera are needed 
only for i and z bands. To create these correction images, 
we first create sky fiat templates. This requires a pro- 
cess of stacking all the detrended images in a band-CCD 
combination after first fiagging all pixels contaminated 
by source fiux. Source contaminated pixels are deter- 
mined by applying a simple threshold above background 
with a variable grow radius so that all neighboring pix- 
els of a pixel determined to contain source fiux are also 
masked. Modal sky values are then calculated for each 
image using pixels that are not fiagged for any reason 
(object pixels, hot column, saturated, interpolated, etc). 
The reduced images are then scaled to a common modal 
sky value, median combined and then rescaled to a unit 
modal value. 

This science sky image then contains a combination of 
any illumination and fringe signatures that are common 
to the input images. To create the illumination correc- 
tion we adaptively smooth the science sky images with 
a kernel that is large in the center and grows smaller 
near the edges. This effectively averages out the effects 
of any fringing, leaving an illumination correction image 
behind. The fringe correction is then produced by first 
differencing the science sky image and the illumination 
correction image, leaving behind an image of the small 
scale structure (i.e. fringe signature) that is common to 
all the science images. This fringe correction image is 
then scaled by the model value of the science fiat image 
to produce a fractional fringe correction image. 

The illumination correction image is applied like a 
fiat field correction to all previously corrected images, 
thereby removing any trends that are introduced by the 
differences in illumination of the dome fiats and the fiat 
sky. The fringe correction is applied by first scaling the 
correction image by the modal value of the sky in the sci- 
ence image and then subtracting it. The results of these 
two corrections are visually very impressive. The fringe 
effects in i and z band are nicely removed in almost all 
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Fig. 5. — Distortion map produced by SCAMP for one Mosaic2 
exposure consisting of 8 images. TPV distortion model was used. 
The Mosaic2 distortions were modelled for each CCD by expressing 
distortions along the RA and DEC direction each with a third order 
polynomial in CCD x and y. 

cases. We have found some problem images where the 
fringe correction leaves clearly visible fringe signatures, 
and these are cases where only a few frames in z or z 
were taken on a particular night, and the use of images 
from neighboring nights to create the corrections was not 
adequate. 

We expect that the residual scatter we measure could 
be further reduced using a star flat technique to bet- 
ter characterize the non- uniformities in the pupil ghost. 
Nevertheless, the delivered data quality from our current 
flattening prescription produces data that meet our data 
quality requirements. We note that the same prescrip- 
tion has been used previously to meet the data qual- 
ity requirements of the SuperMACHO experiment in the 
processing of Mosaic2 data. 

At the end of this series of image detrending steps 
which includes overscan, bias, flat-field, pixel-scale, il- 
lumination and fringe corrections, the pipeline creates 
eight images (one for each CCD) for every science expo- 
sure. These single epoch image FITS files are called red 
images, and they contain 3 extensions: the main image, 
a bad pixel mask (BPM) and an inverse variance weight 
image. The BPM contains a short integer image where 
any unusable pixels have non-zero values (coded accord- 
ing to the source of the problem). The weight-map is 
an inverse variance image map that tracks the noise on 
the pixel scale and where the weight is set to zero for all 
masked pixels. 

3.1.3. Astrometric Calibration 

Besides pointing errors, wide-field imagers exhibit 
instrumental distortions that generally deviate signifi- 
cantly from those of a pure tangential projection. In 
addition, the vertical gradient of atmospheric refractiv- 
ity creates a small image flattening of the order of a few 
hundredths of a percent (corresponding to a few pixels on 
a large mosaic) , with direction and amplitude depending 
on the direction of the pointing. These three c ontribu- 
tions are modeled in the SCAMP fBertinl2006') package 
that we use for astrometric calibration. SCAMP uses 
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Fig. 6. — Median value of the difference in RA and DEC for ob- 
jects in BCS coadd catalogs vs USNOB catalog for every tile in 
arcsecs. The matching is done in a 2" window. The histograms 
are peaked at ~ 0.0104^^ and 0.0084^^ in A-RA and A-DEC respec- 
tively. The rms of the histograms in RA and DEC is about 0.047 
and 0.045 arcsecs. Note that th e intrinsic a ccuracv of the USNOB 
catalog is about 0.2 arcsecs ( ^Monet et al.|2003J . 

the TPV distortion model^^, which maps detector coor- 
dinates to true tangent plane coordinates using a poly- 
nomial expansion. 

SCAMP is normally meant to be run on a large set of 
S Extractor catalogs extracted from overlapping ex- 
posures together with a reference catalog, in order to 
derive a global solution. However, since our pipeline op- 
erates on an image- by-image basis, we proceed in two 
steps: we first run SCAMP once on a small subset of 
catalogs extracted from BCS mosaic images to derive an 
accurate polynomial model of the distortions where the 
distortions in RA/DEC tangential plane are expressed as 
a third degree polynomial function of the CCD x/y posi- 
tion. This Mosaic2 distortion map modeled using a third 
order polynomial per CCD for a BCS exposure is shown 
in Fig. [5] The astrometric solution computed in this 
first step of calibration is based on a set of overlapping 
catalogs from dithered exposures which provides tighter 
constraints on non-linear distortion terms (than catalogs 
taken individually). Using this model, we create a dis- 
tortion catalog that encodes the fixed distortion pattern 
of the detector. We then run SCAMP on catalogs from 
each individual exposure (i.e. the union of the catalogs 
from each of the eight single epoch detrended images), 
allowing only linear terms (two for small position off- 
sets and four for the linear distortion matrix) describing 
the whole focal plane to vary from exposure to exposure. 
The solutions for the World Coordinate System (WCS) 
including the TPV model parameters are then inserted 

-"^^ currently under review for inclusion 
in t he registry of FITS conventions; see 
|http: //fits. gsfc. nasa.gov/registry/tpvwcs .html 
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TABLE 1 
Sextractor Detection Parameters 



Parameter 



Values 



DETECT_TYPE 

DETECT_MINAREA 

DETECT_THRESH 

ANALYSIS_THRESH 

FILTER 

FILTER_NAME 

DEBLEND_NTHRESH 

DEBLEND_MINCONT 

CLEAN 

CLEAN_PARAM 
BACKPHOTO THICK 



CCD 
5 

1.5 
1.5 

Y 

gauss_3 . 0_3x3 . conv 
32 
0.005 
Y 
1.0 
24.0 



back into image headers. This approach capitahzes on 
the expected constancy of the instrumental distortions 
over time. 

We use the USNO-Bl (M onet et aI]|2QQ3| catalog 
as the astrometric reference, tor astrometric refine- 
ment, the cataloging is done using SExtractor, and 
using WINdowed barycenters to estimate the positions 
of sources. 

The astrometric accuracy is quite good, as can be 
demonstrated with the BCS coadds. First, the accuracy 
is at the level of a fraction of a PSF or else significant PSF 
distortions would appear in the coadds, and this is not 
the case. Second, we can measure the absolute accuracy 
relative to the calibrating catalog USNOB by probing 
for systematic offsets in RA or DEC between our object 
catalogs and those from the calibration source. Fig. [6] 
shows the distribution of median offsets within all the 
coadd tiles for both RA and DEC. The mean of the his- 
tograms is 0.0104'' in RA and 0.0084'' in DEC, and the 
corresponding rms scatter is 47 milli-arcsec and 45 milli- 
arcsec, respectively. The USNOB catalog itself has an 
absolute accuracy wi th characteristi c uncertainty of 200 
milli-arcseconds (Monet et al.'"2003), which then clearly 
dominates the astrometric uncertainty of our final cata- 
logs. 

3.1.4. Single- Epoch Cataloging 

To catalog all objects from single-epoch images we 
run SExtractor using PSF modeling and model-fitting 
photometry. A PSF model is derived for each CCD image 
using the PSFEx package ( |Bertin|20lT| . PSF variations 
within the each CCD are modeled as a A^^^ degree poly- 
nomial expansion in CCD coordinates. For our applica- 
tion we adopt a 26 x 26 pixel kernel and follow variations 
to 3^^ order. An example of variation of the Full- Width 
at Half Maximum (FWHM) of the PSF model across a 
single-epoch image is shown in Fig. [7| The FWHM varies 
at the 10% level across this CCD due to both instrumen- 
tal and integrated atmospheric effects. 

A new version of SExtractor (version 2.14.2) uses 
this PSF model to carry out PSF corrected model fitting 
photometry over each image. The code proceeds by fit- 
ting a PSF model and a galaxy model to every source 
in the image. The two-dimensional modeling uses a 
weighted that captures the goodness of fit between the 
observed fiux distribution and the model and iterates to 
a minimum. The resulting model parameters are stored 
and "asymptotic" magnitude estimates are extracted by 
integrating over these models. This code has been exten- 



PSF FWHM 



Fig. 7. — Variation of the PSF model FWHM for g-hsmd across 
a single-epoch image from the BCS night 20061030. Variations 
across the roughly 10^ x 20' image are at the 10% level. 



sively tested within the DESDM program on simulated 
images, but the BCS data provide the first large scale real 
world test. For the BCS application we adopt a Sersic 
profile galaxy model that has an ellipticity and orienta- 
tion. This model fitting is computationally intensive and 
slows the "lightning-fast" SEXTRACTOR down to a rate 
on the order of 10 objects/s on a single core. The SEx- 
tractor config file detection parameters are shown in 
Table m 

The advantages of model fitting photometry on single 
epoch images that have not been remapped are manifold. 
First, pixel to pixel noise correlations are not present in 
the data and do not have to be corrected for in estimat- 
ing measurement uncertainties. Second, unbiased PSF 
and galaxy model fitting photometry is available across 
the image, allowing one to go beyond an approximate 
aperture correction to aperture magnitudes often used 
to extract galaxy and stellar photometry. Third, there 
are morphological parameters that can be extracted af- 
ter directly accounting for the local PSF, which allows 
for improvements in star-galaxy classification and the ex- 
traction of PSF corrected galaxy shear. A more detailed 
description of these new SExtractor capabilities along 
with the results from an extensive testing program within 
DESDM will appear elsewhere (Bertin et al, in prepara- 
tion). 

3.1.5. Remapping 

From the WCS parameters which are computed for ev- 
ery reduced image, one can approximate the footprint of 
the CCD on the sky using frame boundaries in Right As- 
cension and Declination. For the BCS survey we have 
a predefined grid of 36' x 36' tangent plane tiles cov- 
ering the observed fields. Based on this, for every red 
image which is astrometrically calibrated, we determine 
which tiles it overlaps. We then use SWarp (Bertin e t al.| 
2002) to produce background-subtracted remapped im- 
ages that conform to sections of these tangent plane tiles. 
A particular red image can be remapped to up to four dif- 
ferent remap images in this process. Pixels are resampled 
using Lancz6s-3 interpolation. 

Remapping also produces a pixel weightmap and we 
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also remap the bad pixel map (using nearest neighbor 
remapping). In this process of remapping, zero weight 
pixels in the reduced images generically impact multiple 
pixels in the remap image given the size of the inter- 
polation kernel. These remaps are then stored for later 
photometric calibration and coaddition. This on-the-fly 
remapping need not be done, because at a later stage of 
coaddition one could in principle return to the red im- 
ages, but given the PSF homogenization we do prior to 
coaddition we have found it convenient to do the remap- 
ping as we are processing the nightly datasets. 

3.1.6. Nightly Photometric Calibration 

Our initial strategy for photometric calibration in- 
volved traditional photometric calibration using the stan- 
dard fields observed on photometric nights along with the 
image overlaps to create a common zeropoint across all 
our tiles. In fact, within DESDM we have developed a 
so-ca lled Photometric Standards Module (j Tucker et al. 
20071 ) (PSM) that we use to fit for nightly photometric 
solutions, and then we apply those solutions to all science 
images and associated catalogs from that night. For BCS 
this involves determining the zeropoints of all images on 
photometric nights through calibration to identified non- 
variable Standard stars from SDSS Stripe 82 field (Smith 
et al.||2002D. 



Tms procedure was used for processing and calibra- 
tion of the BCS data processed in Spring 2008. But 
closer analysis of these data showed that we were not 
able to control photometric zeropoints to the required 
level to allow for cluster photometric redshifts over the 
full survey area. We therefore abandoned this method 
for BCS in favor of relative photometric calibration us- 
ing common stars in overlapping red images followed by 
absolute calibratio n usin g the stellar locus (described in 
more detail in Sec. 3.2.5). One problem we faced is that 



so-called photometric nights exhibited non-photometric 
behavior in the standard field observations. There was 
no reliable photometric monitor camera at CTIO during 
our survey, and so observers simply used the time hon- 
ored tradition of watching for clouds to make the call 
on a night being photometric. Because of our strategy 
for standard star observations (beginning, middle, end 
of night), even those nights that exhibit consistent pho- 
tometric solutions need not have been photometric over 
the full nights. Therefore, we felt it safer to assume that 
no night was truly photometric and to calibrate the data 
using an entirely different approach. 

The results from the PSM module for those nights ex- 
hibiting good photometric solutions are still useful. They 
have allowed us to monitor changes in the detectors and 
to measure the color terms in transforming our photom- 
etry onto the SDSS system. We provide a brief descrip- 
tion of this procedure, although no science results in this 
paper are based on PSM related direct photometric cal- 
ibration. We expect to apply this method for absolute 
photometric calibration of DES data where we will in- 
deed have an IR photometric monitoring camera on the 
mountain. The PSM solves for the following equation : 

^inst — ^std = an-\-bnX{stdColor — stdcoloro)-\-kX (2) 

where is the photometric zeropoint for all 8 CCDs, bn 
is the color term, stdColor is the fiducial color around 
which we define our standard solutions, which is g — r for 



g and r bands, and r — i for i and z bands, stdcoloro is a 
constant equal to ^ — r = 0.53 for ^ and r bands and r — 
i = 0.09 for i and z bands, k is the first-order extinction 
coefficient and X is the airmass. The PSM module solves 
for a^, bn and k for each photometric night. Using these 
values for the PSM a, b and /c, one can also estimate the 
expected zeropoint for every exposure. We calculate it 
as follows 



ZP = -a + 2.5 log (exptime) - kX 



(3) 



We applied the PSM on about 30 BCS nights which were 
classified as photometric. We also checked for trends in 
variation of color terms as a function of CCD number. 
Only the i-band color term shows some variation, and 
this approximate constancy of color terms greatly sim- 
plifies the coaddition of the data, because we don't have 
to track which CCDs have contributed to each pixel on 
the sky. The color terms we have used for photometric 
calibration are -0.1221, -0.0123, -0.1907, and 0.0226 in 
griz respectively. We also examine the band dependent 
extinction coefficient (k) calculated using data from the 
photometric nights. For the ensemble of about 30 pho- 
tometric solutions in each band, we find the median griz 
extinction coefficients at CTIO over the life of the survey 
to be 0.181, 0.104, 0.087 and 0.067 mag./airmass respec- 
tively. 

This completes the description of all the steps of the 
nightly processing or single-epoch processing that we do 
for BCS. 

3.2. Coaddition 

Once we have data processed for all of the BCS nights, 
we then combine data within common locations on the 
sky to build deeper images that we call coadds. This 
process is called co-addition and is complicated because 
it involves combining data taken in widely separated 
times and under very different observing conditions. Co- 
addition processing is done on a tile by tile basis. We 
describe our approach below. 

3.2.1. Relative Photometric Calibration 

During single-epoch processing we extract instrumen- 
tal magnitudes. To produce science ready catalogs, we 
must calculate the zeropoint for every image and re- 
calibrate the magnitudes. The photometric calibration 
is done in two steps. The first step is a relative zeropoint 
calibration that uses the same object in overlapping ex- 
posures, and the second is an absolute calibration using 
the stellar locus. 

The relative calibration is done tile by tile rather than 
simultaneously across the full survey. We use two dif- 
ferent pieces of information to calculate the relative ze- 
ropoints. The primary constraint comes from the aver- 
age magnitude differences from pairs of red images with 
overlapping stars. The stars are selected based on the 
SEXT RACTO R ffags and spreadjnodel (discussed later 
in Sec. 3.2.4) values. In cases where there aren't enough 



overlapping stars, we use the average CCD to CCD ze 
ropoint differences derived from photometric nights. In 
previous versions of the reduction we also used dir ect ze- 
ropoints derived from photometric nights (see Sec. 3.1.6) 
and relative sky brightnesses on pairs of CCDs. As previ- 
ously mentioned, the direct photometric zeropoint infor- 
mation is contaminated at some level. The sky brightness 
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constraints also seem to be problematic for BCS, because 
only g and r band data were taken on dark nights with 
no moon present which can introduce a gradient across 
the camera. To avoid a degradation of the calibration we 
used neither the sky brightness constraints nor the direct 
photometric zeropoints. 

We determine the zeropoints for all images in a tile 
by doing a least squares solution using the inputs de- 
scribed above. For this least squares solution there are 
N input images, each with an unknown zeropoint in the 
vector z. We arbitrarily fix the zeropoint for one image 
and calibrate the remaining images relative to it. We 
have M different constraints in the constraint vector c. 
The matrix AisNxM and denotes the images involved 
for each constraint. The resulting system of equations is 
described hy Az = c where we use singular value decom- 
position to solve for the vector z. This gives the relative 
zeropoints needed to coadd the data for a particular tile. 

3.2.2. PSF Homogenization 

Combining images with variable seeing generically 
leads to a PSF that varies discontinuously over the coad- 
ded image. This affects star galaxy separation and con- 
tributes to variation across the image in the complete- 
ness at a given photometric depth. The PSF accuracy 
could be quite poor in regions where there are abrupt 
changes to the PSF which would translate into biases 
in the photometry that would be difficult to track. The 
main steps involved in the process of PSF homogeniza- 
tion include: (1) modeling the PSF using PSFEx for all 
remap images contributing to a coadd tile, (2) choosing 
the parameters of the target PSF, (3) using PSFEx to 
generate the homogenization kernel, and (4) carrying out 
the convolution to homogenize all the remap images to a 
common PSF. 

To reduce PSF variation we processed our images to 
bring them to a common PSF within an image and from 
image to image within a coadd tile. To do this we ap- 
ply position dependent convolution kernels that are de- 
termined using power spectrum weighting functions that 
adjust the relative contributions of large scale and small 
scale power within an image in such a way as to bring the 
PSFs within and among the image samples into agree- 
ment. The target PSF is defined to be a circular Moffat 
function with the FWHM set to be the same as the me- 
dian value of all input PSFs. 
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where Yi are the elements of a polynomial basis in x — y. 
The target PSF is defined to be a circular Moffat func- 
tion with the FWHM set to be the median FWHM of 
the input images. We imposed a cut on input image 
PSF FWHM < 1.6 arcsec. This selects only images with 
relatively good seeing. Images from each band are ho- 
mogenized separately. The FWHM of the target PSF for 
all BCS tiles is shown in bottom panel of Fig. [S] 

Another price of homogenization is that noise is cor- 
related on the scale of the PSF. While the noise is al- 
ready correlated to some degree through the remapping 
interpolation kernel, PSF homogenization characteristi- 
cally affects larger angular scales than does the remap- 
ping kernel. This leads to biases in photometric and 
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Fig. 8. — FWHM of single epoch images using PSFex (top panel) 
along with the target PSF FWHM used for homogenizing the coadd 
images for the full BCS survey (bottom panel). The peak values 
of target PSFs are about 1" for g and r bands, 0.9^^ for z-band and 
0.8^^ for z-band respectively. 



morphological uncertainties, and can also affect initial 
object detection process in S Extractor. To address 
this within DES we account for the noise correlations 
on two critical scales by producing two different weight 
maps. The first weight map is used to track the pixel 
scale noise, and the second weight map is used to correct 
for the correlated noise on the scale of the PSF. The pixel 
scale weight map is used by S Extractor in determining 
photometric and morphological uncertainties. The PSF 
scale weight map is used by S Extractor in the detec- 
tion process. Extensive tests within DESDM have shown 
this approach to be adequate to produce unbiased pho- 
tometric and morphological uncertainties and to enable 
unbiased detection of objects within coadds built from 
homogenized images. These results will be presented in 
detail elsewhere. For the BCS processing we used only a 
single pixel-scale weight map, tuned to return the correct 
measurement uncertainties within S Extractor. 

3.2.3. Stacking Single Epoch Images 

We use SWarp to combine the PSF-homogenized im- 
ages to build the coadd tile. Inputs include the rela- 
tive fiux scales derived from the calibration described in 



Sec. 3.2.1). We combine the homogenized remap images 
using the associated weight maps and bpm for each im- 
age. The values of the fiux-scaled, resampled pixels for 
each image are then median combined to create the out- 
put image. This allows us to be more robust to transient 
features such as cosmic rays in the i and z bands where 
there are three overlapping images. Also, objects with 
saturated pixels in all single epoch images will contain 
pixels that are marked as saturated in the coadd im- 
ages as well. This ensures accurate fiagging of objects 
with untrustworthy photometry during the coadd cata- 
loging stage. The resulting output coadd image's size is 
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Fig. 9. — Variation of the PSF model FWHM for a ^-band coadd 
image for the coadd tile BCS0516-5441. Because of the homoge- 
nization process the variations are at the level of 1% across the 36^ 
image. 



8192 X 8192 pixels or approximately 0.6 x 0.6 degrees. 

Fig. [9] shows a map of the FWHM as a function of po- 
sition over one homogenized coadd image. Variations are 
at the level of ~1% over the coadd, as compared to the 
~10% variations that are typical for Mosaic2 across a sin- 
gle CCD (see Fig.[7|. The constancy of PSF as a function 
of a position ensures that the PSF model can be modeled 
accurately and that the PSF corrected model fitting pho- 
tometry is unbiased. The PSF homogenization process 
also circularizes the PSF. Fig [lO] shows the distribution 
of ellipticities for the Mosaic2 single epoch images and 
coadded images (color coded by band). The single epoch 
ellipticity varies up to 0.1 with a modal value around 
0.02. By contrast, the ellipticity distribution of the BCS 
coadds is peaked at a fraction of a percent with a median 
value of 0.001. 

3.2.4. Cataloging of Coadded images 

To catalog the objects from coadded images, we run 
SExTRACTOR in dual-image mode with a common de- 
tection image across all bands. For BCS we use the i- 
band image as the detection image, because it has three 
overlapping images so the cosmic ray removal is good, 
and it is by design the deepest of the bands. We then 
run SEXTRACTOR using model-fitting photometry using 
this detection image and coadded image in each band. 
This ensures that a common set of objects are cataloged 
across all bands. In both single epoch and coaddition cat- 
aloging, the detection criterion was that a minimum of 5 
adjacent pixels had to have flux levels about 1.5<j above 
background noise. The full SExtractor detection pa- 
rameters used for both coadded and single-epoch images 
are shown in Tab.[T] In all we catalog about 800 columns 
across four bands. However for the public data release, 
we have released 60 columns from SExtractor per ob- 
ject. This full list in can be found in Tab.[3j Most of the 
parameters are described in the SExtractor manual 
online. There are a few additional parameters which are 
not yet released in the public version of SExtractor. 
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Fig. 10. — Mean Ellipticity calculated by PSFEx for single-epoch 
images (top panel) and for PSF homogenized coadds (bottom 
panel), color coded by band. Ellipticity is defined as (a — b)/(a -\- b) 
where a and b refer to semi-major and semi-minor axis respectively. 
For single-epoch images the median ellipticty for griz bands is 
0.0342, 0.0326, 0.0374 and 0.04 respectively. For coadds, typical 
values are around 0.004, 0.0024, 0.0026, 0.0033. 



These include model magnitudes and a new star-galaxy 
classifier called spreadjnodel, which is a normalized sim- 
plified linear discriminant between the best fitting local 
PSF model (0) and a slightly more extended model (G) 
made from the same PSF convolved with a circular ex- 
ponential disk model with scale — length = FWHM/ 16 
(where FWHM is the full-width-half maximum of the 
PSF model). It can be defined by the following equation 



spreadjnodel = 



(5) 



where x is the image vector centered on the source. The 
distributio n of s preadjnodel for BCS catalogs is dis- 
cussed in 3.2.7 More details of spreadjnodel will be 



described elsewhere (Bertin et al in preparation). 

3.2.5. Absolute Photometric Calibration 

Once all objects from the coadd are cataloged (in in- 
strumental magnitudes), we proceed to obtain the ab- 
solute photometric calibration using Stellar Locus Re- 
gression (High et al. 2009). The principle behind this is 
that the regularity of the stellar main sequence leads to a 
pre-determined line in color-color space called the stellar 
locus. This stellar locus is observed to be invariant over 
the sky, at least for fields that lie outside the galactic 
plane. The constancy of the stellar locus has been used 
as a cross check of the p hotometric calibration within the 
SDSS survey (Ivezic et a l. 2007). 

Absolute photometric calibration is done after the end 
of coaddition. We select star-like objects using a cut on 
the SExtractor spreadjnodel parameter and magni- 
tude error. We then match the observed stars to 2MASS 
stars from the NOMAD catalog, which is a combination 
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Fig. 11. — The stellar locus in three different color-color spaces 
for the BCS tile BCS0510-5043. The blue hue shows the expected 
distribution derived from studies of a large ensemble of stars within 
the SDSS and 2MASS surveys. Red points show model magnitudes 
of stars from the BCS catalogs of this tile. The stellar locus distri- 
butions allow us to calibrate the absolute photometry and to assess 
the quality of the photometry for each tile. 



of USNOB and 2MASS datasets and which have JHK 
magnitudes ( Skrutskie et al.|2QQ6 ). Color offsets are var- 
ied until the observed locus matches the known locus. 
Because the 2MASS magnitudes are calibrated with a ze- 
ropoint accuracy at the ~2% level, one can bootstrap the 
calibration to the other bands. The known locus is de- 
rived using the high q uality ^^superclean" SDSS-2MASS 
matched catalog from Covey et al.| (2007). It consists of 
^300,000 high quality stars with data m ugrizJHK . A 
median locus is calculated for each possible color combi- 
nation in bins of g — i. 

The fit is done in two stages. First, a three paramater 
fit is done to the g — r — i — z colors. Another 
fit is done using g — r and r — J where the shift in the 
g — r color is fixed from the first fit. The first fit provides 
an accurate calibration of the colors and the second fit 
fixes the absolute scale. The fit is done this way because 
only a fraction of the stars that overlap with 2MASS are 
saturated in all bands. We perform the stellar locus cali- 
bration for model and 3 arcsec aperture magnitudes, sep- 
arately. The model magnitude calibration is then used to 
calibrate the other magnitudes in the catalog (except for 
the 3 arcsec magnitudes). The calibration of the 3 arcsec 
aperture magnitudes determines a PSF dependent aper- 
ture correction for the mag_aper_3 magnitudes only. We 
have found these small aperture magnitudes to provide 
higher signal to noise colors for faint galaxies in compar- 
ison to mag_model and mag_auto. 

An example of the stellar locus fits for one coadd tile is 
shown in Fig. [TT] Red points are the observed colors of 
the stars, and the blue line is the median SDSS-2MASS 
locus. The orthogonal scatter about the stellar locus for 
all 3 color combinations for BCS tiles is shown in Fig 
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The rms orthogonal scatter about the stellar locus in 
[g — r^r — i), {r — i^i — z) and {g — r, r — j) is 0.059, 0.061 
and 0.075 respectively. Given the scatter and the number 
of stars available for calibration we can determine the 
zeropoints in our bands with sub-percent accuracy. 

3.2.6. Testing Stellar Locus Calibration in SDSS 

To validate our photometric calibration algorithms, we 
applied exactly the same procedure to the full SDSS- 
2MASS catalog in ( [Covey et al.||2007D . This catalog in- 
cludes noiser objects than the catalog we used to derive 
the median stellar locus. We selected 4 areas (between 
RA of 120° and 350°) and divided each into 1° x 1° 
patches. We match the objects to obtain 2MASS mag- 
nitudes and then apply the same calibration procedure 
as we did for the BCS catalogs. The rms scatter dis- 
tributions for all the three color combinations can be 
found in Fig. [12] (bottom panel). The corresponding 
scatter for SD^ in (^ — r, r — z), (r — — z), and 
(^ — r, r — j) is about 0.041, 0.035 and 0.05 respectively 
and is about 1.5 times smaller than for the BCS catalogs. 
This is clear evidence for higher scatter in our stellar pho- 
tometry as compared to the SDSS photometry. Assum- 
ing this additional source of scatter adds in quadrature 
with the SDSS observed scatter, we estimate the extra 
noise in BCS color combinations compared to SDSS is 
0.039, 0.054 and 0.048 in (^ - r, r - i), {r-i,i-z), and 
{g — r^r — j) respectively. Because these noise sources are 
getting contributions from each color, we can estimate 
that the noise floors are S{g — r) ~ 0.027, S{r — i) ~ 0.038, 
5{i — z) ^ 0.038. These then imply noise floors in the 
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Fig. 12. — Stellar locus scatter (above) for three color combina- 
tions for all tiles in BCS survey (top panel) and the same for SDSS 
(bottom panel). Typical BCS scatter is in the 5% to 8% range, and 
offsets after calibration are characteristically 1% or less. Typical 
scatter and offsets in the SDSS dataset are smaller than in the BCS 
survey, reflecting the tighter requirements on photometric quality 
in SDSS. 



stellar photometry within griz bands of approximately 
1.9%, 2.3%, 2.7% and 2.7%, respectively. This is then 
in good agreement with the typical repeatability scatter 
seen in these bands (see Fig. 14) when one considers that 
g and r bands each have two overlapping exposures and 
i and z band each have three. 



3.2.7. Star- Galaxy Classification 

Our current catalogs contain two star-galaxy classifica- 
tion parameters provided by SExtractor : class_star 
which has been extensively studied and spreadjnodel, 
which has been newly developed as part of the DESDM 
development program. In order to test their performance 
and range of magnitudes up to which these measures can 
be reliably used, we plot the behavior of these two classi- 
fiers in i-band as a function of mag-model in Fig. 
class_star lies in the range from to 1. At brigEt 
magnitudes one can see two sequences in class_star 
for galaxies and stars near and 1 respectively. The 
two sequences begin merging as bright as i = 20 and 
are significantly merged beyond i = 22. As described 
in Sec. 3.2.4| spreadjnodel uses the local PSF model 
to quantify tne differences between PSF-like objects and 
resolved objects. In the spreadjnodel panel it's clear 
that there is a strong stellar sequence around the value 
0.0, and that galaxies exhibit more positive values. The 
narrow stellar sequence and the broad galaxy sequence 
begin merging at z = 22 in the BCS, but there is signif- 
icant separation in the two distributions of points down 
to i = 23. Along with spreadjnodel comes a measure- 
ment uncertainty, and so it is, for example, possible to 
define a sample of objects that lie off the stellar sequence 
in a statically significant way. For the BCS data, a good 
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Fig. 13. — Plots of spreadjnodel (top panel) and class_star 
(bottom panel) as a function of i-band magnitudes for the full BCS 
catalog. Note that both measurements exhibit separate sequences 
for stars and galaxies, and that as one moves to fainter magnitudes 
these sequences merge. This is simply due to low signal to noise 
objects not containing enough morphological information for a re- 
liable classification. However, note also that the new spreadjnodel 
retains good capability of separating galaxies from stars to fainter 
magnitudes than class_star. 

cut to separate stars would be spreadjnodel < 0.003. 
Detailed studies of this new classification tool have been 
carried out within the DESDM project and will be car- 
ried out elsewhere. 

3.2.8. Quality Control and Science Ready Catalogs 

During the processing within the DESDM system a va- 
riety of quality checks are carried out. These include, for 
example, thresholding checks on the fraction of flagged 
pixels within an image and the ^ind number of stars 
used in the astrometric fit of each exposure. In addition, 
the system is set up to report on the similarity between 
correction images (bias, flat, ilium, and fringe) against 
stored templates that have been fully vetted. During the 
BCS processing this last facility was not used. 

Our experience has been that problems at any level of 
processing are most likely to show up in the stages of rela- 
tive and absolute photometric calibration. Therefore, for 
the BCS processing done here we capture a range of pho- 
tometric quality tests including the number of stars used 
in the stellar locus calibration and the rms scatter about 



the true stellar locus of the calibrated data (see Fig. 11). 
In addition, we examine the photometric repeatabilityTor 
common objects within overlapping images contributing 
to each tile. In Fig. [M] we show an example for the g- 
band in tile BCS0549-5043. This shows the magnitude 
difference between pairs of overlapping objects versus the 
average magnitude. The scatter here includes both sta- 
tistical and systematic contributions, and the envelope 
of scatter grows toward faint magnitudes, as expected. 
Outlier rejection is done on the point distribution, and 
all 3a outliers are flltered out and colored red. In the top 
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Fig. 14. — Repeatability plots for single-epoch images for BCS 
tile BCS0549-5043 in gf-band. The repeatability is used to test the 
quality of the photometry in each band and tile. The top panel 
shows the mean magnitude difference between different single- 
epoch images which cover same region of sky, binned as a function 
of magnitude along with statistical and systematic errors. The bot- 
tom panel shows an un-binned representation of the same. Char- 
acteristic scatter on the bright end (i.e. the systematic floor) is 2% 
to 3% for g and r and 3% to 4% for i and z. 



panel we plot the mean and rms as well as the outlier 
fraction of these repeatability distributions as a function 
of magnitude. The mean and rms numbers are listed 
in milli-mags. Also, the statistical uncertainties of the 
model magnitudes are used to estimate the systematic 
magnitude error contribution to the rms. On the bright 
end where the statistical noise is very small, the system- 
atic contribution to the rms is close to the total, which is 
10 millimags in this case. As one moves toward the faint 
end the statistical contribution increases and the esti- 
mated systematic contribution plays only a small role in 
explaining the scatter. This is just how we expect the 
photometry to behave. 

Repeatability plots indicate systematic contributions 
to the photometric errors at the 10 to 20 milli-mag levels 
for typical g and r band tiles. For i and z band tiles 
the systematic noise is closer to 30 to 40 milli-mags. For 
all the BCS tiles we have examined these repeatability 
and stellar locus plots to probe whether the scatter is 
in acceptable ranges. In cases where tiles didn't meet 
these quality control tests we worked on the relative and 
absolute photometric calibration to improve the data. 
In addition to these photometry tests we examined the 
sky distribution of cataloged objects within each tile. In 
cases where large numbers of faint "junk" objects were 
found we attempted to remove them by adjusting the 
cataloging. At present all our tiles meet these quality 
tests except for a handful of tiles that are marked in red 
in Fig. [l] This includes 4 tiles in the 5 hr field and 6 tiles 
in the 23 hr field, corresponding to ~4% of the 80 deg^ 



region. Ideally, we would reimage these regions to obtain 
better data. 

For every BCS night, the detrending pipeline creates 
three main types of science image files which we denote 
as raw^ red and remap. The coadd pipeline produces four 
coadd images per tile for each of the four bands. Once we 
have calibrated coadd catalogs for all the processed tiles 
we run a post-processing program to remove duplicate 
objects near edges of the tiles. This is necessary because 
there is a 2 arcmin overlap between neighboring tiles. 
The program selects sources that appear in neighboring 
tiles that lie within in 0.9 arcsec radius and for each pair 
it keeps the object that lies farthest from the edge of 
its tile. In this way a single, science ready catalog is 
prepared for each field. The 23 hr field catalog contains 
1,877,088 objects, and the 5 hr field contains 2,952,282 
objects with i model magnitude < 23.5. In the next 
section we review additional tests of the data quality. 
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Fig. 15. — Histogram of lOcr magnitude limits for all BCS tiles 
using mag.auto errors in all four bands. The median depth values 
for all BCS tiles are 23.3, 23.4, 23 and 21.3 in griz, respectively. 
The corresponding lOcr point source depths are 23.9, 24.0, 23.6 and 
22.1. 



3.3. Survey Depth 

We estimate the lOa photometric depths for galaxies 
using SEXTRACTOR mag.auto errors. This is obtained 
by doing a linear fit to the relationship between the mag- 
nitude and the log of the inverse magnitude error. As a 
cross-check we also estimated the depths using informa- 
tion in the weight maps, and the results were comparable. 

The distributions of depth for each band over the full 
survey is shown in Fig. [15] The median magnitude 
depths for griz bands are 23.3, 23.4, 23.0 and 21.3 re- 
spectively. These numbers are shallower than the depths 
we estimated using the NOAO exposure time calculator 
during the survey planning; those depths were 24.0, 23.9, 
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23.6, and 22.3. Our originally proposed depths assume a 
2.2 arcsec diameter aperture, whereas galaxies near the 
10a detection threshold are typically larger in our im- 
ages. We examine the depths of 2^' aperture photometry 
and find that the median depths are 24.1, 24.1, 23.5 and 
22.2 in griz^ respectively. These are within 0.2 mag of 
our naive estimates, explaining the bulk of the differ- 
ence. In addition, we know that during our survey often 
the conditions were not photometric, and this could in- 
troduce another 0.1 to 0.2 mag offset. Another reason 
for the difference is that the calibrated observed magni- 
tudes also include a correction for galactic extinction and 
reddening, whereas the estimated depths did not have 
extinction corrections included. 

Corresponding lOa point source depths are extracted 
using model fitting mag-psf uncertainties. The results 
in bands griz are 23.9, 24.0, 23.6 and 22.1, respectively. 
These are in better agreement with the small aperture 
photometry depths we used to estimate the exposure 
times for the survey. 

Another way of probing the depth of the survey is 
to look at the number counts of sources as a function 
of magnitude. Fig. [16] contains the logN-logS from the 
combined 5 hr and 23nr fields using mag.auto. No star- 
galaxy separation is carried out, because near the detec- 
tion limit there is not enough morphological information 
to reliably classify. The magnitudes of the turnover in 
the counts corresponds to 24.15, 23.55, 23.25 and 22.35. 
These turnover magnitudes mark the onset of significant 
incompleteness in the catalogs. Estimates of the depth 
of the 50% and 90 % completenesss limi t for a subset of 
the tiles appear in Zenteno et al. (2011), but we do not 
apply that analysis to the whole survey. 
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Fig. 16. — Number counts of BCS objects for all four bands in 
the BCS field using mag -auto. The turnover magnitudes are 24.15, 
23.55, 23.25 and 22.35 in griz, respectively. The corresponding 
median mag-auto 10a depths in griz are 23.0, 23.4, 23.0 and 21.3 
respectively. 



Finally, we probe for spatial variations in photometry 
by examining the distribution of sources above certain 
flux cuts over the two survey regions. The distributions 
of all sources at i < 22.5 in both the 5 hr and the 23 hr 
fields are shown in Fig.[T7[ Objects are excluded for tiles 
that did not pass our quality tests, and this produces 4 
black squares in the 5 hr field, and 6 black squares in 
the 23 hr field. The general uniformity of this object 
density distribution is an indication that the absolute 
photometric calibration is reasonably consistent across 
the fields. In the 5 hr field it is clear that one tile in 
the lower left does not reach the depth i = 22.5 reached 
by the other tiles. This defect disappears if we examine 
the density distribution roughly 1 magnitude brighter, 
indicating this is a depth issue and not a photometric 
calibration problem. For the 23 hr field there is a small 
black rectangular notch in the upper right of the field 
with an associated dark path. Within this tile we have 
verified that too few of our z-band exposures met the 
seeing requirements, and that has led to an uncovered 
region (the notch) as well as the shadow of lower object 
density to the right. Again, this is a depth issue rather 
than a photometric issue. There is another shadowed tile 
visible in the lower right part of the 23 hr field, and this 
is also a depth issue. 

In Fig.[T8]we show similar object density plots for stars 
and galaxies for the 5 hr field. The stars and galaxies 
are chosen based on spreadjnodel cut at 0.003, where 
all objects with values greater than this threshold were 
considered galaxies and the rest were considered stars. A 
catalog depth cut at z < 22.5 was imposed. This is shown 
in Fig. [is] The stellar distribution is quite uniform across 
this fiel37 indicating that spreadjnodel performance is 
quite robust to variations of PSF across a survey. Note 
that the shallow tile in the lower left portion of the survey 
exhibits edge effects, which we believe are associated with 
the reduced depth of this tile relative to the others. In the 
lower panel is the galaxy distribution. The same shallow 
tile shows up in the lower left portion. In addition, it is 
clear that the galaxy density is varying as a function of 
position as expected for the large scale structure of the 
Universe. We are quite happy with this performance. We 
have explored the same plots in the 23 hr field, and the 
results are similar. Moreover, we have explored these 
plots created using star_class as the classifier. The 
spatial distribution is highly inhomogeneous, indicating 
that class_star cannot be used to reliably separate stars 
and galaxies in a uniform manner across a large survey. 

3.4. BCS Data release 

We are publicly releasing the BCS catalogs, images 
and the photo-z training fields to the astrophysical com- 
munity. All public BCS data products can be down- 
loaded from http://www.usin.uni-inuenchen.de/BCS, 
The BCS catalogs are divided into ascii files for the 5 
hour and 23 hour fields. Separate catalogs are available 
for the tiles that passed our quality analysis and for the 
tiles that did not. Each catalog contains 63 columns 
which are described in Table |3] We are also making 
available the coadded images for the BCS survey at the 
same site. These images are available in a PSF homog- 
enized form (used for the cataloging) and in the non- 
homogenized form. As in the case of the catalogs, we 
split the tar files by field and by whether the tiles passed 
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Fig. 17. — Distribution of sources in 5 hr (top) and 23 hr fields 
(bottom) from the combined catalogs after a mag-auto magnitude 
cut (i < 22.5). The gaps show tiles which were not included in 
the release due to data quality problems. Some other tiles have 
only partial coverage or do not push to the depth of the magnitude 
cut with good completeness. A logarithmic scale (zscale option 
in ds9) is used. The uniformity of the source distribution is a 
demonstration of the photometric uniformity across the survey. 



our quality tests or did not. These tarballs contain FITS 
tile compressed images, which reduces the volume by a 
factor of ~5 relative to the uncompressed coadds. 

4. PHOTOMETRIC REDSHIFTS 

Initial tests of data quality are undertaken by obtain- 
ing photometric redshifts for BCS objects using an arti- 
ficial neural network. Neural networks have been used to 
determine accurate photomet ric redshift s in past optical 
surveys (Collister et al., 2QQ7| [Oyaizu et al.,.2QQ8bJ . We 




Fig. 18. — Distribution of stars (top) and galaxies (bottom) in 
BCS 5 hr field with mag-auto i < 22.5 based on spreadjnodel 
cut of 0.003. The stars look uniformly distributed and you can see 
traces of large-scale structure in galaxy density plot. A logarithmic 
scale (zscale option in ds9) is used. We have explored similar plots 
with class_star, and these contain very large inhomogeneities in 
the stellar and galaxy distribution, indicating that spreadjnodel 
offers significant advantages over class_star in the classification 
of objects in large surveys. 



use ANN2;, a feed-forward multi-layer perceptron network 
designed for finding photometric redshifts (Colliste r fc| 
'Lahav||2QQ4 ) . The network is composed of a series ot m- 
puts, several layers of nodes, and one or more outputs. 
Each node is made of a function that takes its input as 
a weighted output of each of the previous layer's nodes. 
The weights are tuned by training the network on a rep- 
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Fig. 19. — Comparison between calibrated model magnitudes for 
stars from four BCS standard tiles (after stacking them together) 
and SDSS magnitudes after color and extinction corrections for g- 
band. The stars are chosen by requiring that class_star > 0.8 in 
all four bands and also SExtractor flag < 5. The histograms are 
normalized to unity. The peak offset between BCS model magni- 
tudes and SDSS is —0.06 in r, i and 0.02 in z bands. 



resent ative dataset with known outputs. The optimal set 
of weights are those that minimize a cost function, which 
reflects the difference between a known output value and 
the network's predicted value. 

The training process can result in a set of weights that 
are over-fit to a particular training set. Furthermore, a 
given training process can converge to a local minimum of 
the cost function instead of the true minimum. In ANN2;, 
the first issue is overcome by finding the set of weights 
that minimizes the cost function on a separate validation 
set rather than on the training set itself. The second is 
avoided by training a committee of several networks with 
randomized initial weights. The mean weights from each 
committee are used in the final network. 

We train our neural network on 5,820 objects with 
known redshifts. It is run with eight input parameters: 
four magnitudes griz] three colors g — r — and i — z; 
and a concentration index. Mag.auto magnitudes are 
used for individual filters, mag-aperS magnitudes are 
used to determine colors, and the i-band spreadjnodel 
is used for the concentrat ion i ndex. Following the guide- 
lines of [Firth et al] ( |2003D and |Collister fc Lahav| |2004D , 
we use a minimally sutticient network architecture and 
committee size in hopes of achieving the highest quality 
results. We find this to be a committee of 8 neural net- 
works that each have an architecture of 8:16:16:1 (eight 
inputs, two hidden layers of 16 nodes each, and one out- 
put). We denote photometric and spectroscopic redshifts 



as Zphot and 

^phot 



-'spec5 



respectively and have Az represent 



We compare our photometry with SDSS data by look- 
ing at spectroscopic calibration tiles which overlap with 
SDSS data and which contain significant numbers of 
spectroscopic redshifts. As explained in Sec. |4.2[ these 
spectroscopic redshifts are then used for training our neu- 
ral networks to obtained photometric redshifts. To do a 
comparison with SDSS catalogs, we applied color and 
extinction corrections to SDSS catalogs from these tiles. 
The fields which we consider for this purposes are from 
CNOC and DEEP fields centered at RA, DEC values of 
(02 hr 25 min, 7°), (02 hr 29 min, 35°), (23 hr 27 min, 
8°), and (23 hr 29 min, 12°). 

We then do an object by object comparison of colors 
and magnitudes of all stars from SDSS versus those from 
BCS catalogs in these tiles. The SDSS magnitudes for 
objects which overlap BCS tiles go up to 23.4 in g and 
21.6 in r, z, and z. We consider an object to be matched 
if it spatially overlaps to within 2" . Since the number 
of objects in each tile which overlap with SDSS is small, 
we combine results from all tiles into one plot for each 
magnitude or color as necessary. The magnitude com- 
parison for all four bands (using mag -model) is shown in 
Fig.[T9j The peak offset between BCS model magnitudes 
and^DSS magnitudes is approximately —0.06 in g^r and 
i bands and about 0.02 in z bands, while the median 
offset is —0.0562 in ^, r and i and 0.0087 in z. 

We also do a color comparison using the same cuts 
for these tiles between BCS and SDSS colors using 
mag-aperS (magnitude within a 3 arcsec aperture), be- 
cause colors are determined using this magnitude in 
photo-z estimation (Fig.[2Q|). The peak offset in colors in 
g—r^r—i^ and i—z is about —0.01, —0.03, and —0.02 mag- 
nitudes respectively. The median offset is about —0.01 
for ^ — r and i — z and about —0.05 for r — i. The rms 
scatter about the median is 0.052, 0.061 and 0.081 for 
g — r — and i — z^ respectively. 

4.2. Photometric Redshift Calibration 

We obtain our training dataset by dedicating nine of 
the survey pointings to fields overlapping spectroscopic 
surveys: CDFS, CFRS, two CN0C2 fields, SSA 22, three 
DEEP2 fields, and VVDS. Objects from these fields share 
their photometric depth and reduction pipeline with the 
BCS data as well as have known spectroscopic redshifts. 
Although this training set is not representa tive of the 
survey in sky position, Abdalla et al. (20111) show that 
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4.1. Photometry Crosschecks with SDSS 



limiting a neural network trainmg set to small patches 
of sky does not result in biased redshifts for large sur- 
veys. The key issue is having uniformity of photometry 
between the training and application fields. 

Only objects that have reliable redshifts and photo- 
metric parameters are used to train the neural network. 
Objects with an i-band magnitude > 22.5 or an i-band 
error > 0.1 are removed from the training set. Objects 
that are unresolved in one or more bands or that have 
a SExtractor flag greater than 2 are removed as well. 
Similar cuts are made based on spectroscopic redshift er- 
rors, however the nature of the cut varies by catalog. The 
DEEP2, CN0C2, and CFRS catalogs provide redshift er- 
rors for each measurement. Objects from these fields are 
removed if their spectroscopic redshift errors are greater 
than 0.01. The ACES catalog (providing coverage of the 
CDFS field) and the VVDS catalog assign a confidence 
parameter for each object. In this case we only include 
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Fig. 20. — Difference in (g — r), (r — i), and (i — z) colors between 
stars from BCS tiles and SDSS using mag-aperS. All cuts are same 
as in Fig. |19|and histograms are normalized to unity. 



TABLE 2 
Photo-z Training Fields 



Survey 


RA 


Dec 


Redshifts 


ACES^ 


03:32 


-27:48 


2846 


CFRS^^ 


22:17 


00:91 


65 


CNOC2^ 


02:25 


00:07 


318 


CNOC2^ 


02:26 


00:43 


164 


SSA 22^ 


01:40 


00:01 


818 


DEEP2^ 


02:29 


00:35 


226 


DEEP2^ 


23:27 


00:08 


414 


DEEP2^ 


23:29 


00:12 


600 


VVDS^ 


14:00 


05:00 


329 



Cooper et al.|(|201ll 


Lilly et al. 


(|l995b|a| 


Yee et al.|i 


200011 



^ Cowie et al | (|l994i 



I Davis et al.|(|2 007); "New man et al.| ( |2012| ) 
Le Fevre et a l. (2005 200^ 

objects with a confidence of 3 or 4 (see respective sur- 
veys for definitions). Both primary and secondary tar- 
gets from the VVDS survey are included. Additional cuts 
were experimented with but found to produce more out- 
liers, a larger sigma, or to reduce the size of the training 
set too much. 

The final training set contains 5,820 objects. Table [2] 
breaks down the number of training objects that pass 
the filter criteria from each pointing. Figure [2l] further 
breaks down these objects by redshift bin. The pointings 
combine to provide a consistent distribution of redshifts 
from < 2: < 1.1. 

We have released the matched catalogs of spectroscopic 
redshifts along with information from BCS catalogs for 
these fields. This would enable others to develop their 
own photometric redshift estimates using these data. 
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Fig. 21. — Redshift distribution of 5,820 objects from the cali- 
bration fields used to train ANNz. The redshift distribution is color 
coded by source. 



We evaluate the performance of ANN2; on our data by 
randomly selecting half of the objects from the training 
set to train ANN2; while the other half remains for test- 
ing. One sixth of the objects from the training half are 
removed to form the validation set (see above). The re- 
sult provides 2,910 objects with both photometric and 
spectroscopic redshifts. 

We measure the photometric redsh ift performance us - 
ing three metrics. The first, following Ilbert et al. (2006), 



is the normalized median absolute deviation 

\Az\ 



(^Az/{l+z) 



1.48 X median 



spec 



This metric is better suited for our data than the stan- 
dard deviation as it is less affected by catastrophic out- 
liers. The second is the fraction of catastrophic outliers 
T] defined as the percentage of objects that satisfy 



\Az\ 



l + Zs 



> 0.15. 



(6) 



The third metric is the net bias in redshift, averaged over 
all N objects and defined as 



^bic 



1 ^ 



TV 



Our training set performs as a^. 



7(1+^) 



= 0.061 with 



T] = 7.49%. Over the entire range of redshifts there is lit- 
tle net bias: Zbias = 0.0005. These statistics, particularly 
the fraction of catastrophic outliers, can be improved by 
culling objects based on their photometric redshift error. 
ANN2; provides redshift errors that are derived from the 
errors of the input parameters, however there are sev- 
eral ot her methods of determining photometric redshift 
errors. Oyaizu et al.| ([2008a) evaluate how well various 
methods improve 2;phot statistics. They show that culling 
objects based on redshift errors derived from magnitude 
errors are competitive with other methods at reducing 
the redshift scatter and catastrophic outlier fraction. 

We analyze our photometric performance after culling 
our data of objects with Zphot error > 0.13 based on er- 
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Fig. 22. — Top panel: Two-dimensional histogram of Zphot vs. 
Zspec for training set objects that have 2;phot error < 0.13. Bin 
sizes are O.OlSz x O.OlSz. Red bins count catastrophic outhers as 
defined by Equation [6] Blue bins count all other objects. 5 objects 
with Zspec or Zp]^ot > 1-5 are not displayed. Bottom panel: The 
same training set data are shown with each point representing a 
bin of 50 objects. 



rors provided by ANN2;. The performance of the culled 
data improves to (Jaz/{i-\-z) = 0.054 and rj = 4.93%. 
However, the 2;bias increases slightly to 0.0022. While 
the bias increases, it is still negligible. Figure 22 demon- 
strates the performance of ANN2; in determining redshifts. 
For objects within the range 0.3 < z^pec ^ 0.9, our photo- 
metric and spectroscopic redshifts match with little bias. 
For objects with redshifts below 0.3, there is a positive 
bias and for objects with redshifts beyond ^spec ~ 0.9 
there is a negative bias. 



4.3. Application to the full BCS Catalog 

The 5,820 objects from the training set are used to 
train a committee of 8 ANN2; networks, each with an ar- 
chitecture 8:16:16:1. This committee is used to deter- 
mine redshifts and errors for every object in the BCS 
catalog. These are included in columns 61 and 62 of 
the data release. Because we found a negligible net bias 



when testing our calibration set, we do not perform a 
bias correction to redshifts of the BCS catalog. 

Many of the objects of the BCS catalog lie outside of 
the parameter space of data used to train ANN2;. While 
Collister fc Lahavj ( |2004 ) have demonstrated success us- 
ing ANN2: to determme redshifts of galaxies outside the 
parameter space used to train the network, this was done 
using a set of galaxies with a very uniform distribution 
of spectral types. For the generic distribution of galaxies 
provided in the BCS catalog, neural networks are unreli- 
able in predicting redshifts outside the trained parameter 
space. Therefore we indicate whether an object lies in- 
side or outside of the parameter space of the training set 
with a flag (column 63). A value of 1 means the ob- 
ject is within the parameter space of the training set and 
the redshift is reliable. A value of means the object 
lays outside the parameter space and the redshift is un- 
reliable. The flag is based only on the magnitude and 
magnitude error cuts that were made on the training set 
(i.e. i < 22.5, z-error < 0.1, resolved in all bands). It 
is not based on the SExtractor flag, the star-galaxy 
separation criteria, or on photometric redshift errors. 

We were able to obtain photo-z's for about 1,955,400 
objects from the BCS catalog with i < 22. From these, 
there are ~204,600 objects in the catalog that pass the 
star-galaxy separation criteria in all bands and lie within 
the training set parameter space. The redshift distribu- 
tion of these BCS objects in different magnitude ranges 
is shown in Figure |23j The peak redshift is around 
^phot = 0-4 for 20 < z < 22. Out of these, there are 
about 200 objects with 2:phot > 1-0. 

Many of the objects that do not pass star-galaxy sepa- 
ration are stars. Since ANN2; was trained with only galax- 
ies, stellar objects lie outside the parameter space and 
therefore do necessarily not get assigned a correct red- 
shift of 2;phot = 0. In fact, only a handful of objects 
in the entire catalog are assigned a %hot close to zero. 
We investigate the performance of the photo-z's when 
training a network with both stars and galaxies. Using 
the same inputs and network architecture as above but 
including '^1,000 stars in the training set, ANN2; was suc- 
cessful in assigning stars a redshift below 0.1 only 70% 
of the time. However, redshift assignment of galaxies 
was not adversely affected. Only 4 out of approximately 
3,400 galaxies were assigned a redshift less than 0.015. 
The fraction of catastrophic outliers as well as o-/\z/{i-\-z) 
were not significantly affected either, so long as stars 
are not included in the statistics. While the results of 
training ANN2; with stars are not sufficient to use for the 
entire BCS catalog, th ese preliminary result s show some 
promise. Furthermore, Collister et al. (2007) have shown 



better results when training an ANN2; network specifically 
for star-galaxy separation. 

5. CONCLUSIONS AND DISCUSSION 

In this paper we present an overview of the Blanco Cos- 
mology Survey (BCS), an ~80 deg^ optical photometric 
survey in griz bands carried out with the Mosaic2 im- 
ager on Blanco 4m telescope between 2005 — 2008. We 
discuss the observing strategy within the context of our 
scientific goals, and we present basic observing character- 
istics at CTIO such as the sky brightness and delivered 
image quality. 

We provide a detailed description of the data process- 
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Fig. 23. — Distribution of photometric redshifts of the galaxies 
that he within the ANNz training set parameter space and have 
^phot error < 0.05(1 + 2;phot) and pass the star-galaxy separation 
test and in different i-band model magnitude ranges. 



ing, calibration and quality control, which we have car- 
ried out using a development version of the Dark En- 
ergy Survey Data Management System. The processing 
steps in going from raw exposures to science ready cat- 
alogs include image detrending and astrometric calibra- 
tion; this processing is run independently on every night 
of observations. This is followed by image co- addition, 
which combines data from the same region of the sky into 
deeper coadd images. 

The processing of real data from the Blanco telescope 
provides a real world stress test of the DESDM sys- 
tem. Many novel algorithmic features, which will be 
used to process upcoming DES data, were tested on BCS 
data. These include PSF homogenization, cataloging us- 
ing PSF corrected, model-fitting photometry, object clas- 
sification using the new spreadjnodel, absolute photo- 
metric calibration using the stellar locus and a variety of 
quality control tests. 

We present the characteristics of the dataset, includ- 
ing the median estimated lOa galaxy photometry depth 
in the coadds for bands griz^ which are 23.3, 23.4, 23.0, 
and 21.3, respectively. The corresponding point source 
10a depths in griz are 23.9, 24.0, 23.6, and 22.1, re- 
spectively. We measure the systematic noise fioor in 
our photometry using photometric repeatability in sin- 
gle epoch images and comparisons of the stellar locus 
scatter from BCS and SDSS. Both results indicate a 
noise fioor at the ~ 1.9% level in ^, ~ 2.2% in r, and 
~ 2.7% in i and z bands. This noise fioor does not im- 
pact the core galaxy cluster science for which the BCS 
was designed. We expect that with an improved char- 
acterization of the illumination correction using the star 
fiat technique demonst rated in the CFHT Legacy Sur- 
vey fRegnaul t et al.|2Q09| ) it would be possible to reduce 
this noise fioor lurtner, but ffiven that the current fioor 



is adequate for our science needs we have not included 
these corrections in our BCS processing. 

Our absolute photometric calibration is obtained using 
the stellar locus and including the 2MASS J band pho- 
tometry. We can calibrate our zeropoints at better than 
~1% (statistical) to the stellar locus, and so our overall 
photometric uniformity is driven by the ^2% acc uracy 
of the 2MASS survey (e,g, [Skrutskie et al.||2006D . We 
show that our photometric zeropoint calibration is quite 
uniform across the survey by showing star and galaxy 
counts across the survey. We also demonstrate that with 
spreadjnodel it is possible to carry out uniform star- 
galaxy separation even across a large extragalactic sur- 
vey. 

As an additional data quality test, we present photo- 
metric redshifts derived from a neural network trained on 
a sample of objects with spectroscopic redshifts that we 
targeted during the BCS survey. The performance of our 
four band griz photometric redshifts are evaluated based 
on analysis of a calibration set of over 5,000 galaxies with 
measured spectroscopic redshifts. We find good perfor- 
mance with a characteristic scatter of crA2/(i+2) ~ 

0.054 

and an outlier fraction of 7^ = 4.93%. Finally, we pro- 
vided a summary of the output data products from our 
co-added images and catalogs along with information on 
how to download them. 

Finally, the BCS data have been used for a range of 
scientific pursuits, which we briefiy summarize and ref- 
erence here to allow the reader to seek additional in- 
formation as needed. Within the SPT survey, the first 
four SZE selected clusters were optic ally confirmed with 
redshift estimates using BCS data (Staniszewski et al. 



|2009|) and detailed studies of galaxy population s usin j 
these clusters were reported in Zenteno e t al.] (2011). 
The total number of SPT cluster candidates wit h si^nal- 



to-noise rat io > 4.5 in BCS footprint is 15 (Reichardt 



et al. 2012|) and among these 10 have been confirmed 



with the BCS data and the r emaining 5 have redshift 
lower limits between 1 and 1.5 fSong et al.||2012 ). These 
clusters and their BCS derived redshifts hav e figured 
prominently in SPT publications to date ([S taniszewski] 
eraI]|2QQ9l |Vanderlinde et al.|[2QTQl [High et al.H2010i 



ct di. ^yjyjc/ V ciiiu.c± iiiiU-C cu cii.|| ±±i^ii cii.h^vjxvj} 

Andersson et al.||2011| IW ilham son et al.|2011[ |Reichardt| 
et al.||2012'; 'Song et al.( |2012| ). The BCS data enabled 



the serendipitous discovery of a strong lensing arc of a 
galaxy at z = 0.9057 b y a massive galaxy cluste r at a 
redshift of z = 0.3838 ( |Buckley-Geer et al.|[20lT] ). Ad- 
ditional automated searches lor strong lensing arcs have 
also been carried out, and further analysis of BCS data 
for weak lensing is in progress. 

A sample of about 105 galaxy clusters was found us- 
ing first thr ee seasons of BC S data u sing a n independent 
processing ([Menanteau et a l. 2009, 2010b), and the BCS 
data were also used for optica l confirmation of ACT clus- 
ters ( Menanteau et aL]|2010a ). Other stu dies include es- 



timat es ol weak lensing cluster masses ( Mclnnes et al 
2009|) and a search for Q SO candidates using r-band 



data (Jimenez et al. 2009). 

We used the BCS data to measure photometric red- 
shifts of abo ut 46 X-ray selected clusters in the XMM- 
BCS survey ("Suhada et al. 2012). This X-ray selected 
sample is currently being used in combination with SPT 
data to explore the low mass cluster population and its 
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SZE properties (Liu et al., in prep). In addition, these 
BCS data are also being used in the analysis of the larger 
XMM-XXL survey in the 23 hr field (Pierre, private com- 
munication). 

The BCS data continue to provide an important 
dataset for SPT. Recently, the data were used to trace 
the galaxy populations and were correlated against th e 
SPT CMB-lensing maps ( |van Engelen et al.| [20121 ), 
demonstrating correlation s signiticant at th e 4a to 5cr 
level in both BCS fields (Bleem et al.||2012D. The BCS 



data will provide a valuable optical dataset for combina- 
tion with a 100 deg^ Spitzer survey over the same region 
(Stanford, private communication), a 100 deg^ Herschel 
survey (Carlstrom, private communication), and they 
will overlap one of the deep mm-wave fields being tar- 
geted by SPT-pol (Carlstrom, private communication) 
until the Dark Energy Survey data are available. 



ing us spectroscopic redshifts for objects from the SSA 
22 field. The Munich group acknowledges the support of 
the Excellence Cluster Universe and from the program 
TR33: The Dark Universe, both of which are funded by 
the Deutsche Forschungs Gemeinschaft. We acknowledge 
support from the National Science Foundation (NSF) 
through grants NSF AST 05-07688, NSF AST 07-08539, 
NSF AST 07-15036, and NSF AST 08-13534. We ac- 
knowledge the support of the University of Illinois where 
this project was begun. This paper includes data gath- 
ered with the Blanco 4-meter telescope, located at the 
Cerro Tololo Inter- American Observatory in Chile, which 
is part of the U.S. National Optical Astronomy Observa- 
tory, which is operated by the Association of Universities 
for Research in Astronomy (AURA), under contract with 
the NSF. 

Facilities: Blanco (MOSAIC) 



We would like to acknowledge Len Cowie for provid- 



REFERENCES 



Abdalla, F. B., Banerji, M., Lahav, O., & Rashkov, V. 2011, 

MNRAS, 1460 
Abell, G. O. 1958, ApJS, 3, 211 

Abell, G. O., Corwin, Jr., H. G., & Olowin, R. P. 1989, ApJS, 70, 
1 

Adami, C., et al. 2010, A&A, 509, A81+ 
Andersson, K., et al. 2011, ApJ, 738, 48 

Bert in, E. 2006, in Astronomical Society of the Pacific Conference 
Series, Vol. 351, Astronomical Data Analysis Software and 
Systems XV, ed. C. Gabriel, C. Arviset, D. Ponz, & S. Enrique, 
112-+ 

Bertin, E. 2011, in Astronomical Society of the Pacific Conference 
Series, Vol. 442, Astronomical Data Analysis Software and 
Systems XX, ed. I. N. Evans, A. Accomazzi, D. J. Mink, & 
A. H. Rots, 435-+ 
Bertin, E., Mellier, Y., Radovich, M., Missonnier, G., Didelon, P., 
&; Morin, B. 2002, in Astronomical Society of the Pacific 
Conference Series, Vol. 281, Astronomical Data Analysis 
Software and Systems XI, ed. D. A. Bohlender, D. Durand, Sz 
T. H. Handley, 228-+ 
Birkinshaw, M. 1999, Physics Reports, 310, 97 
Biviano, A. 2008, ArXiv e-prints, 0811.3535 
Bleem, L. E., et al. 2012, ArXiv e-prints, 1203.4808 
Bruzual, G., &; Chariot, S. 2003, MNRAS, 344, 1000 
Buckley-Geer, E. J., et al. 2011, ApJ, 742, 48 
Butcher, H., &; Oemler, Jr., A. 1984, ApJ, 285, 426 
Carlstrom, J. E., Holder, G. P., &; Reese, E. D. 2002, ARA&A, 
40, 643 

Cioni, M.-R. L., et al. 2011, A&A, 527, A116 
Collister, A., et al. 2007, MNRAS, 375, 68 
Collister, A. A., &; Lahav, O. 2004, PASP, 116, 345 
Cooper, M. C, et al. 2011, ArXiv e-prints, 1112.0312 
Covey, K. R., et al. 2007, A J, 134, 2398 

Cowie, L. L., Gardner, J. P., Hu, E. M., Songaila, A., Hodapp, 

K.-W., & Wainscoat, R. J. 1994, ApJ, 434, 114 
Cunha, C, Huterer, D., & Dore, O. 2010, Phys. Rev. D, 82, 

023004 

Dalai, N., Dore, O., Huterer, D., & Shirokov, A. 2008, 

Phys. Rev. D, 77, 123514 
Davis, M., et al. 2007, ApJ, 660, LI 

de Jong, J. T. A., Verdoes Kleijn, G. A., Kuijken, K. H., 
Valentijn, E. A., KiDS, & consortiums, A.-W. 2012, ArXiv 
e-prints, 1206.1254 

Dressier, A. 1980, ApJ, 236, 351 

Firth, A. E., Lahav, O., & Somerville, R. S. 2003, MNRAS, 339, 
1195 



Fowler, J. W. 2004, in Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, Vol. 5498, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference 
Series, ed. C. M. Bradford, P. A. R. Ade, J. E. Aguirre, 
J. J. Bock, M. Dragovan, L. Duband, L. Earle, J. Glenn, 
H. Matsuhara, B. J. Naylor, H. T. Nguyen, M. Yun, &; 
J. Zmuidzinas, 1-10 

Gal, R. R., Lopes, P. A. A., de Carvalho, R. R., Kohl-Moreira, 
J. L., Capelato, H. V., &; Djorgovski, S. G. 2009, AJ, 137, 2981 

Gilbank, D. G., Gladders, M. D., Yee, H. K. C, & Hsieh, B. C. 
2011, AJ, 141, 94 

Gladders, M. D., Yee, H. K. C, Majumdar, S., Barrientos, L. F., 
Hoekstra, H., Hall, P. B., & Infante, L. 2007, ApJ, 655, 128 

Gonzalez, A. H., Zaritsky, D., Dalcanton, J. J., &; Nelson, A. 
2001, ApJS, 137, 117 

Haiman, Z., Mohr, J. J., & Holder, G. P. 2001, ApJ, 553, 545 

Hansen, S. M., Sheldon, E. S., Wechsler, R. H., & Koester, B. P. 
2009, ApJ, 699, 1333 

High, F. W., et al. 2010, ApJ, 723, 1736 

High, F. W., Stubbs, C. W., Rest, A., Stalder, B., & Challis, P. 

2009, AJ, 138, 110 
Holder, G., Haiman, Z., & Mohr, J. J. 2001, ApJ, 560, Llll 
Ichiki, K., &; Takada, M. 2011, ArXiv e-prints, 1108.4688 
Ilbert, O., et al. 2006, A&A, 457, 841 
Ivezic, Z., et al. 2007, AJ, 134, 973 

Jimenez, R., Spergel, D. N., Niemack, M. D., Menanteau, F., 

Hughes, J. p., Verde, L., &; Kosowsky, A. 2009, ApJS, 181, 439 
Koester, B. P., et al. 2007, ApJ, 660, 239 
Le Fevre, O., et al. 2005, A&A, 439, 845 
. 2004, A&A, 428, 1043 

Lilly, S. J., Hammer, F., Le Fevre, O., &; Crampton, D. 1995a, 
ApJ, 455, 75 

Lilly, S. J., Le Fevre, O., Crampton, D., Hammer, F., & Tresse, 

L. 1995b, ApJ, 455, 50 
Lin, Y., Mohr, J. J., Gonzalez, A. H., & Stanford, S. A. 2006, 

ApJ, 650, L99 

Lin, Y., Mohr, J. J., &; Stanford, S. A. 2003, ApJ, 591, 749 
Lin, Y.-T., & Mohr, J. J. 2004, ApJ, 617, 879 
Majumdar, S., & Mohr, J. J. 2003, ApJ, 585, 603 
. 2004, ApJ, 613, 41 

Mclnnes, R. N., Menanteau, F., Heavens, A. F., Hughes, J. P., 
Jimenez, R., Massey, R., Simon, P., & Taylor, A. 2009, 
MNRAS, 399, L84 

Menanteau, F., et al. 2010a, ApJ, 723, 1523 

Menanteau, F., & Hughes, J. P. 2009, ApJ, 694, L136 

Menanteau, F., et al. 2010b, ApJS, 191, 340 

. 2009, ApJ, 698, 1221 



22 



Desai et al. 



Mohr, J. J., et al. 2008, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 7016, 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series 

Monet, D. G., et al. 2003, AJ, 125, 984 

Newman, J. A., et al. 2012, ArXiv e-prints, 1203.3192 

Ngeow, C, et al. 2006, in Society of Photo-Optical 

Instrumentation Engineers (SPIE) Conference Series, Vol. 6270, 
Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series 

Oyaizu, H., Lima, M., Cunha, C. E., Lin, H., & Frieman, J. 
2008a, ApJ, 689, 709 

Oyaizu, H., Lima, M., Cunha, C. E., Lin, H., Frieman, J., & 
Sheldon, E. S. 2008b, ApJ, 674, 768 

Perlmutter, S., et al. 1999, ApJ, 517, 565 

Regnault, N., et al. 2009, A&A, 506, 999 

Reichardt, C. L., et al. 2009, ApJ, 694, 1200 

. 2012, ArXiv e-prints, 1203.5775 

Rest, A., et al. 2005, ApJ, 634, 1103 

Rozo, E., et al. 2010, ApJ, 708, 645 

Ruhl, J., et al. 2004, in Proc. SPIE, Vol. 5498, Milhmeter and 
Submillimeter Detectors for Astronomy II, ed. J. Zmuidzinas, 
W. S. Holland, & S. Withington (Belhngham: SPIE Optical 
Engineering Press), 11-29 



Schmidt, B. P., et al. 1998, ApJ, 507, 46 

Skrutskie, M. F., et al. 2006, AJ, 131, 1163 

Smith, J. A., et al. 2002, AJ, 123, 2121 

Soares-Santos, M., et al. 2011, ArXiv e-prints, 1111.4717 

Song, J., A, & B. 2012, In prep 

Staniszewski, Z., et al. 2009, ApJ, 701, 32 

Sunyaev, R. A., &: Zel'dovich, Y. B. 1972, Comments on 

Astrophysics and Space Physics, 4, 173 
Tucker, D. L., et al. 2007, in Astronomical Society of the Pacific 

Conference Series, Vol. 364, The Future of Photometric, 

Spectrophotometric and Polarimetric Standardization, ed. 

C. Sterken, 187-+ 
Suhada, R., et al. 2012, A&A, 537, A39 
van Engelen, A., et al. 2012, ArXiv e-prints, 1202.0546 
Vanderlinde, K., et al. 2010, ApJ, 722, 1180 
Wang, L., &; Steinhardt, P. J. 1998, ApJ, 508, 483 
Williamson, R., et al. 2011, ApJ, 738, 139 
Yee, H. K. C, et al. 2000, ApJS, 129, 475 
Zenteno, A., et al. 2011, ApJ, 734, 3 



APPENDIX 
BCS CATALOG DESCRIPTION 



We created an ASCII catalog files which is obtained from the catalogs of each individual tile and after removing 
duplicates. The description of each column in the BCS catalog is provided in Tab. |3] 
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TABLE 3 
Details of BCS Catalogs 



Column 


Parameter 


Units 


Definition 


1 


"tilsnsiiii© 






Name of tile 


2 


ob j ectid. 






ID from DESDM database table coadd-objGcts 


3 


RA 


de^ 


^ree 


Right Ascension 


4 


DEC 


de^ 


^ree 


Declination 


5 


in a.g jh d e 1 _g 


AB 


mag 


Model magnitude (^) 


6 


mag G rr jnod G 1 _g 


AB 


mag 


Error in model magnitude (^) 


7 


nia.g_a.ut; o _g 


AB 


mag 


Kron magnitude (^) 


8 


magG 3r3r -aut o _g 


AB 


mag 


Error in Kron magnitude 


9 


ma.g,_poj. _g, 


AB 


mag 


PSF magnitude (p) 


10 




AB 


mag 


n.T*TT»T* in P.^Ti TTi n orn 1'1'nH p> ( n\ 

J— /1 1 Ill J. iiidgiiit Line yy J 


11 




AB 


mag 


Po-l-Tr~»c;i n n tti n crni'l"nrip> i n\ 
J. c; ti *-»oiciii iiidgiii t LHJ.C \}J ) 


12 




AB 


mag 


Y^iVvcw in Pptrocjipn m a en Hp ( n\ 

l-JLLyJL 111 J. dil l-/olCtll lllCt^lll u LlvlCy \y J 


13. 


m o cr aTioT*^ cr 


AB 


mag 


A/T PI orn 1 "1" n H p in UTT'CPr* T^pyi'iiyp i n\ 
ividgiiiL LIU.C; 111 «j di *^oc*^ djJCitU-ic \[J ) 
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f la.gs_g 






SEXTRACTOR flag (^) 
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SExTRACTOR star/galaxy separator 
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_ 


Difference in PSF and Sersic magnitude 
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Error in model magnitude (t*) 
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Kron magnitude (r) 


22 
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Error in Kron magnitude (r) 
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PSF magnitude (r) 
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Error in PSF magnitude (r) 
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Petrosian magnitude ('^') 
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Magnitude error in 3-arcsec aperture (r) 
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SExTRACTOR flag (r) 
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Petrosian magnitude (i) 
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Magnitude in 3-arcsec aperture {%) 
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SExTRACTOR flag (l) 
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Model magnitude (2;) 
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Kron magnitude (z) 
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Error in PSF magnitude {z) 
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Petrosian magnitude {z) 
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Photometric redshift 
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z_phot_Grr 






Photometric redshift error 
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Within ANN2; training set parameter space 



Note. — Explanation and contents of catalogs in the BCS survey release. More details on some of 
the parameters can be found in the SExtractor manual. The magnitudes are corrected for galactic 
extinction. 



